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Uvodem 

Mile resitelky, mill resitele, 

jiz polovina 32. rocniku je za nami, coz znamena, ze ti nejlepsi z vas mohou v nejblizsi dobe 
ocekavat ve svych schrankach pozvanky na jarni soustredeni. 

Avsak ti z vas, ktcri objevili F'YKOS teprve ted ncbo se vam moc ncdarilo, nczoufcjte! Jiz 
resenim teto scrie si muzete zajistit misto na podzimnim soustredeni, ktere bude urcite skvele! 

Ve ctvrte serii se budeme zabyvat zejmena mechanikou, ale take si vyzkousime, jak moc 
dokaze papir stinit zvuk. V serialu si povime neco o Taylorove rozvoji a ulohach, ktere diky 
nemu dokazeme resit. 

Unor je nabity nejruznejsimi akcemi, za zrninku jiste stoji Fyziklani 2019, ktere se uskutecni 
v patek 15. 2. v Top Ilotelu Praha. Pozor, registrace se uzavira jiz 6. 2.! Pokud se vam nechce 
jezdit do Prahy jen „na otocku“, mame pro vas v ramci Fyziklani pripraveny doprovodny 
program, ktery zacina 14. 2. Jedmm Dnem s Fyzikou na MF'F UK v Praze. V sobotu a v nedeli 
po Fyziklani muzete navstivit radu zajimavych mist a pracovist nejen v ramci MFF UK. 

Zaverem bycliom vain radi jmenem vsecli organizatoru poprali stastny vstup do roku 2019, 
at se vam dari nejen ve skole, ale i v osobnim zivote. Budeme se tcsit na videnou! 

vasi Organizatori 



Zadanf IV. serie 


Termm uploadu: 12. 2. 2019 23.59 
Termm odeslani11. 2. 2019 


Uloha IV. 1 ... kostka se vzduchem 3 body 

Mejme dutou kostku s hranou delky a = 20 cm naplnenou vzduchem s teplotou to = 20 °C, coz 
je zaroven teplota okoli kostky. Vzduch uvnitr kostky ochladime na t\ =5 °C. Jaka sila bude 
pusobit na kazdou stenu kostky? Kostka pri ochlazeni vzduchu v ni nemeni svuj objem. Tlak 
v okoli kostky je po = 101,3 kPa. 


Uloha IV.2 ... utrhne se 3 body 

Mame (nehmotny) provazek delky / a na jeho konci kulicku (hmotny bod) s hmotnosti m. 
Vime, ze maximalm tiha, co unese, je sila F = mg, kde g je mistni tihove zrychleni, ale uz 
nic vie. Provazek upevnime a kulicku budeme drzet ve stejne vysce, jako je misto upevneni, ve 
vzdalenosti delky provazku od druheho konce provazku, ale tak, abychom ho nenapinali. Kulicku 
uvolnime a ta se zacne vlivem tihoveho zrychleni pohybovat. Pod jakym uhlem provazku vuci 
svisle rovine se provazek pretrhne? 


Uloha IV. 3 ... levitujici 6 bodu 

Matej ma rad levitujici veci, a tak si poridil nekonecnou nevodivou nabitou vodorovnou rovinu 
s plosnou nabojovou hustotou a. Pote nad ni umistil micek o hmotnosti m nabity nabojem q. 
Vypocitejte, pro jake hodnoty a muze micek vubec nad deskou levitovat. V jake vysce h se pak 
muze vzna^et? Uvazujte konstantni tihove zrychleni g. 
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Uloha IV.4 . .. trampolma 7 bodu 

Dva hmotne body skakaly na trampoline do vysky ho = 2 m. Vc chvili, kdy oba byly v nejnizsim 
moznem miste trajektoric (vychylka y = 160 cm), jeden z nich zahadne zmizel. Do jako nejvyssi 
vysky byl druhy vymrsten? Kruhova trampolma ma obvod o =10 m a pruzi diky N = 42 pruzi- 
nam s tuhosti A: = 1720 N-m 1 . Trampolmu modelujme N pruzinami rozmistenymi rovnomerne 
a spojenymi ve stredu. Hmotnost zmizeleho hmotneho bodu je M = 400 kg. 

Uloha IV.5 . .. frisbec 9 bodu 

Tenky liomogenni disk obilia na vodorovne podlozce po kruznici s polomerem R. Velikost rych- 
losti teziste disku je v. Urcete uhel a mezi rovinou disku a svislym smerem. Treni mezi diskem 
a podlozkou je dostatecne. Polomer disku je radove mens! nez R. 

Uloha IV.P . . . V-l ve vesimru 10 bodu 

Mezihvezdny prostor neni prazdny, nybrz obsahuje nepatrne mnozstvi limoty. Uvazujte jen vo- 
dfk, potrebnou hustotu si vyhledejte. Molila by existovat kosmicka lod, jez by „nasavala“ vodik 
pfed sebou a vyuzivala energii z nej? Jak rychla/velka by musela byt, aby udrzela termojadernou 
fuzi jen z prijatcho vodiku? Jake jine pfekazky realizacc je rozumne uvazovat? 

Uloha IV.E ... papirova izolacni 12 bodu 

Zmefte, jak moc dokaze papir stfnit zvuk. K mereni stacf pouzft napr. mobilni telefon jako 
generator zvuku a mikrofon v pocitaci jako detektor (Audacity). Pouzijte papiry ruznych druhu 
a tvaru. 

Uloha IV.S ... lagrangeovska 10 bodu 

V zavere serialu ste si urcite vsimli Lagrangian a diferencialnu rovnicu, ktore akoby „spadli 
z neba“. To nie je vobec nahoda, velkou castou tejto serialovej ulohy bude tieto dve rovnice 
odvodit. 

1. Ukazte, ze ak mame pohyb castice v rubovolnom centralnom poli, teda v poli, kde potencial 
zavisi len na vzdialenosti, bude sa castica zarucene pohybovat len v rovine. 

Navod: Zostavte Lagrangeove rovnice II. druhu pre tuto situaciu, pouzite pri tom vhodne 
zovseobecnene suranice. Nasledne bez ujmy na vseobecnosti poloztc suradnicu i? = x/2 
a pociatocnu rychlost v srnere tejto suradnice nulovu. Zamyslite sa a vysvetlite, prcco je 
takato volba v poriadku a nestratJme pri nej ziadne riesenie. 

2. Zostavte Lagrangian pre hmotny bod pohybujuci sa v rovine v centralnom poli. Mali by 
ste dostat ten isty, ako je uvedeny v zavere serialu. Pre tento Lagrangian nasledne najdite 
vsetky intergaly pohybu a pomocou nich najdite diferencialnu rovnicu prveho radu pre 
premennu r. Pre vasu kontrolu, mala by vam vyjst rovnako ako na konci serialu. 

3. Zamyslite sa, ako urcif uhlovu vzdialenost medzi dvoma bodmi na sfere, ak mate zadane 
ich sfericke suradnice. Ukazte to napnklad pre hviezdy Betelgeuze a Sirius, ktorych surad¬ 
nice si najdite. Pomocka: Tato uloha sa da jednoducho vyriesit aj bez znalosti sferickej 
trigonoinetrie. 
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Kategorie ctvrtych rocm'ku 

jmeno 

Student Pilny 

1 . Martin Vanek 

2 . Katerina Rosicka 

3. Ivan Iludak 

4. Jin Blaha 
5.-6. Vojtech Ulman 
5.-6. Leonardo Wimrner 

7. Marie Grunova 

8 . Jindrich Jelinek 

9. Matej Mosko 
10.—11. Michal Juza 
10.—11. Filip Konafik 

12 . Kristian Matustik 

13. Vojtech Klimes 

14. YUQING XU ARIA 

15. Jan Obofil 

16. Tomas Drobil 

17. Jindrich Dusek 

18. Tadeas Wilczek 

19 . Marek Jankola 

20 . Jakub Ruzicka 

21 . Tereza Prokopova 

22 . Katerina Char'vdtovd 

23. Andrcj Rendek 

24. Karel Balej 

25. Miroslav Horsky 

26. Dominik Majkus 

27. Bibiana Hroncova 
28.—30. Bcmadeta Mafikovd 
28.—30. Stepan Tichy 

28.—30. Marie Vanakovd 
31.—32. Vratislav Besta 
31.—32. Marek Talir 
33. Marek Bozon 


skola 1234 5PES III % E 

MFF UK 3 3 6 8 8 10 12 10 60 100 180 


G, Vysoke My to 4 3 6 8 8 - - 10 39 93 127 

G J. Ortena, Kutna Ilora 3 2 6 8 10 3 - 7 39 88 107 

ESS, Lip. Mikulas 3 3 6 - - - 12 - 24 92 94 

G, Uherske Hradiste 3 3 6 8 - - 12 6 38 82 84 

G Jaroslava Seiferta, Praha 2364 145 - 25 ^7 75 

Colegio pH, Tijuca, Brazil 4 3 6 8 6 3 10 7 47 85 75 

G, Moravsky Krumlov 1 2 -- - 10 13 58 74 

G, Olomouc-Hejcin - ____ — 108 65 

G Grosslingova, Bratislava 2 3 5 8 - 2 — - 20 70 59 

G, Benesov 3 1 - - - 2 10 - 16 57 54 

G F. Palackeho, Val. Mez. 3 3 6 8 - - - - 20 93 54 

G, Benesov 313-224 - 15 42 51 

G, Tfebon - - - - 67 40 

Dulles High School, USA 4 2 6 8 4 3 4 7 38 63 38 

Klasicke a spanelske G, Br-- — — - — 74 37 

no 

G, Dacice - ____ _ 7 2 36 

G Jana Keplera, Praha - _____ — 57 34 

G F. Zivneho, Bohumin 3 2 6 ---- - 11 73 32 

G M. Ilattalu, Trstena - ____ — 7 1 27 

G, Nymburk - ____ _ 63 2 4 

G Jura Hronca, Bratislava — 1- ____ 1 #5 23 

G B. Nemcove, HK 3 3 6 ----- 12 79 22 

G, Dubnica n. Vahorn - ____ _ 30 18 

G a SOS, Rokycany -_6 - - - - - 6 100 17 

G, Ceskolipska, Praha -- — - - — 52 16 

G Na Vitezne plani, Praha- ____ _ 13 

G Postova, Kosice — 11 - 11 92 11 

G, Cesky Krumlov - ____ _ 56' 9 

G, Jatecni, Usti nad Labem-— — — - — 4 3 9 

G Boticska, Praha - _____ — 55 9 

G, Olomouc-Hejcin -- — — - — 33 6 

G, Cesky Krumlov - _____ _ 55 9 

G, Karvina - ____ _ 3 


39 
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Kategorie tretfch rocm'ku 



jmeno 

Student Pilny 

skola 1 

MFF UK 3 

2 

3 

3 4 
6 8 

5 

8 

P 

10 

E 

12 

S 

10 

Ill 

60 

% 

100 

E 

180 

1 . 

Viktor Matema 

G Brno, tr. Kpt. Jarose 4 

3 

6 

8 

8 

10 

12 

_ 

51 

89 

142 

2 . 

Martin Schmicd 

G, Jihlava 4 

3 

6 

8 

8 

5 

9 

10 

53 

76 

129 

3. 

Jaroslav Herman 

G Brno, tf. Kpt. Jarose 4 

3 

6 

8 

8 

2 

10 

5 

46 

75 

127 

4. 

Jakub Jobus 

G Jura Ilronca, Bratislava 4 

2 

6 

4 

4 

4 

10 

10 

44 

69 

118 

5. 

Jan Benda 

G, Litomericka, Praha 5 

3 

6 

8 

8 

- 

11 

10 

51 

90 

116 

6 . 

Jaroslav Scheinpjlug 

G Jirovcova, Ceske Budejo- 4 
vice 

3 

5 

3 

4 

5 

— 

7 

31 

68 

103 

7. 

Radka Kfizova 

G J. Heyrovskeho, Praha 3 

3 

5 

2 

6 

- 

7 

- 

26 

78 

102 

8 . 

David Komanek 

G, Spitalska, Praha 3 

3 

6 

6 

- 

- 

12 

- 

30 

76 

100 

9. 

Martin Vavfik 

G, Sumperk 4 

- 

5 

- 

10 

- 

- 

- 

19 

87 

95 

10 . 

Jin Zelenka 

G Z. Wintra, Rakovnik 2 

3 

6 

- 

- 

2 

12 

- 

25 

63 

84 

11 . 

Jakub Strnad 

KIvanovo G Kyjov 3 

3 

6 

8 

- 

6 

4 

0 

30 

55 

82 

12 . 

Jan Divila 

G, Lesni ctvrt, Zlin 








— 

52 

52 

13.-14. 

Ronald Dobos 

G Postova, Kosice 








— 

85 

51 

13.-14. 

Pavla Rudolfova 

G, Videnska, Brno 

- 

6 

- 

- 

- 

- 

- 

6 

91 

51 

15. 

Veronika Hendrychova G, Turnov 3 

2 

2 






7 

64 

49 

16. 

Lukas Hronek 

G, Pisek - 







10 

10 

100 

47 

17. 

Jakub Pravda 

SpMNDaG, Bratislava — 








— 

67 

34 

18.-19. 

Katerina Barotova 

G, Olomouc-Hejcin 








— 

60 

28 

18.-19. 

Matej Holubicka 

G, SOS, SOU a VOS, Hori- 3 

2 

- 

- 

- 

3 

6 

- 

14 

56 

28 

20 . 

Matej Kratky 

ce 

PORG, Praha 








_ 

71 

27 

21 .-22. 

Petra Palkovacsova 

G, Nove Zamky 








— 

38 

23 

21 .-22. 

Vaclav Zvomcek 

G Brno, tf. Kpt. Jarose 








— 

72 

23 

23. 

Josef Polasek 

G Jana Keplera, Praha 








- 

48 

22 

24. 

Rachel Johnson 

Richardson High School, - 
USA 








— 

42 

21 

25. 

Milan Tichavsky 

Slezske G, Opava 








- 

67 

20 

26. 

Krystof Jasensky 

G, Lesni ctvrt, Zlin 








— 

45 

19 

27.-28. 

Minh Khoi Ho 

Hanoi - Amsterdam IIS, - 
Vietnam 








— 

76 

16 

27.-28. 

Filip Zukal 

G, Blansko — 








— 

27 

16 

29. 

Eva Vochozkova 

Biskupske G, Brno 








— 

62 

13 

30. 

Jan Vondra 

G Tyn nad Vltavou 1 








1 

48 

12 

31. 

Dominika Kodlovd 

PORG, Praha 








— 

29 

11 

32. 

Daniel Kratky 

G, Trutnov - 








— 

59 

10 

33. 

Katerina Roupova 

G, Blansko 








- 

23 

9 

34. 

Erika Zitniakova 

Evanjclicke G, Banska Bys- - 
trica 








— 

23 

5 

35.-36. 

Lucia Krajcoviechova 

G Jura Hronca, Bratislava — 








— 

67 

4 

35.-36. 

Martina Pivkova 

Evanjclicke G, Banska Bys- — 

- 

- 

- 

- 

- 

- 

- 

- 

40 

4 


trica 
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Resent' III. serie 


Ulolia III.l ... zlevnenc banany 3 body; prumer 2,84; resilo 49 studentu 

Mikulas v obchode vlozil nekolik bananu do igelitovelio sdcku. Pfed jejich zvazemm ho napadlo, 
zc kdyby pytlik naplnil iru'sto vzduchu heliem, budou banany stat o ncco incnc. Helium Mikulas 
koupil vc sieve zajednu korunu na litr pri standardmm tlaku. Jaka musi byt cena bananu, a by 
se mu lento „podvod (C vyplatil? 

Bonus Naleznete plyn, u ktcrelio sc vyplati plnit jun sacck pfi cene bananu 30 korun na 
kilogram. Nezapomchte citovai zdrojc ceny daneho plynu. 

Nad cirri premyslite vy, hdyz v obchode vazite banany? 

Ked’ vrecusko s bananmi naplmme heliom, bude nan posobit vztlakova sila, ktora je podia 
Arcliimedovho zakona rovna F vz = gVg, pricom sa pocita s hustotou tekutiny g, v ktorej je 
objekt ponoreny, a objemorn ponorenej casti objektu V. Zaroven posobf tiazova sila na helium 
samotne, ktora ho „£aha dole“, preto plati 

A rriBfj = (g ™d - ^?He) Vg , 



kde Amo je zmena vahy bananov, g vz d a gn c su hustoty vzduchu a helia a V je objem vo vrec- 
usku, ktory sme vyplnili heliom. Tento objem vierne vypocitat ako podiel hmotnosti pouziteho 
helia a hustoty helia 

v _ WHe 

gue 

cize po vykrateni g a dosadeni V mame vztah 


A ms 


( gw zd 
\ Quo 


1 nil 


Kedze Mikulas chce, aby sa mu podvod opiatil, cena pouziteho helia musi byt minimalne 
rovna usetrenej cene bananov 

rriHeCue < AtubCb , 

kde Cue , Cb su ceny za kilogram helia, resp. bananov. Po dosadeni za A ms dostavame vztah 
pre cenu bananov 

p v C He 

_ l ' 

0 Ha 

Cenu helia v jednotkach CZK-kg _1 vypocitame z hustoty helia a zadanej ceny helia chc 
v jednotkach CZK-^ _1 ako 


Cue = 


c He 
gHe 


Dosadime do finalneho vztahu pre cenu bananov a dostaneme 


C B > -—- = 900 CZK . 

gvzd Pile 
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Kategorie druhych rocm'ku 



jmeno 

skola 

1 

2 

3 

4 

5 

P 

E 

S 

Ill 

% 

E 


Student Pilny 

MFF UK 

6 

6 

6 

8 

8 

10 

12 

10 

GG 

100 

198 

1 . 

Robert Gemrot 

G Komenskeho, Havirov 

6 

6 

6 

8 

9 

10 

12 

4 

61 

80 

150 

2 . 

Martina Dankova 

Klasicke a spanelske G, Brno 

8 

6 

6 

- 

- 

10 

11 

- 

41 

85 

111 

3. 

Patrik Kasparek 

Katolicke gymnazium Treble 

6 

6 

5 

- 

- 

2 

12 

4 

35 

65 

107 

4. 

Soha Husakova 

G, Ceskolipska, Praha 

8 

4 

6 

2 

- 

9 

6 

- 

35 

11 

100 

5. 

Adam Mendl 

G P. de Coubertina, Tabor 

- 

- 

- 

4 

9 

- 

9 

- 

22 

83 

9G 

6 . 

Lubor Cech 

G Brno, tr. Kpt. Jarose 

6 

6 

5 

2 

0 

4 

1 

- 

24 

54 

90 

7. 

Elena Chocholakova 

G L. Svobodu, Humenne 

6 

6 

5 






17 

86 

85 

8 . 

Marco Souza de Joode 

G Nad Stolou, Praha 

6 

6 

5 

- 

- 

2 

- 

6 

25 

66 

84 

9. 

Jan Klivan 

G, Dacice 

4 

4 

5 

- 

- 

- 

9 

- 

22 

63 

64 

10 . 

Tomas Tkileja 

G L. Svobodu, Humenne 

6 

6 







12 

05 

56 

11 . 

Radek Lacko 

G Komenskeho, Havirov 









— 

65 

40 

12 . 

Daniel Perout 

G Brno, tr. Kpt. Jarose 

- 

6 

5 

- 

- 

2 

- 

- 

13 

62 

33 

13. 

Jarmila Tcrpdkovd 

G L. Svobodu, Humenne 









- 

82 

28 

14. 

Karolina Lclochova 

G, Sternbcrk 









- 

51 

26 

15. 

Jan Cervenan 

G J. Pivecky, Slavicin 









- 

52 

23 

16. 

Sarka Stepankova 

G J. Ressela, Chrudim 









- 

63 

19 

17.-18. 

Matej Dvorak 

G Jana Keplera, Praha 









- 

72 

18 

17.-18. 

Lucia Gintnerova 

G Sv. Frantiska, Zilina 









— 

90 

18 

19. 

Martin Polydcsko 

G Alejova, Kosice 









- 

58 

15 

20 . 

Gabriel Scurka 

SpMNDaG, Bratislava 









- 

53 

10 

21 .- 22 . 

Daniel Czinege 

SPS chemicka, Ostrava 









— 

47 

9 

21 .- 22 . 

Kristyna Chlupacovd 

G J. Ressela, Ghrudim 

0 

4 

5 






9 

50 

9 

23. 

Vojtech Janota 

G, Strakonicc 









— 

47 

7 

24. 

JYantisek Krus 

Masarykovo G, Plzen 









— 

27 

6 

25. 

Laura Ilanouskovd 

G J. Wolkera, Prostejov 

- 

2 

- 

- 

- 

1 

- 

- 

3 

19 

3 


37 
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Poradi resitelu po III. serii 


Kompletm vyslcdky najdotc na http://fykos.cz. 


Kategorie prvm'ch rocm'ku 



jmeno 

skola 1 

2 

3 

4 

5 

P 

E 

S 

Ill 

% 

£ 


Student Pilny 

MFF UK 6 

6 

6 

8 

8 

10 

12 

10 

66 

100 

198 

1 . 

Eva Feldbabelova 

Katolicke gymnazium Treble 0 

6 

6 

6 

2 

- 

- 

7 

27 

77 

128 

2 . 

Jim Kohl 

Biskupske G, Brno 6 

6 

5 

8 

6 

- 

12 

- 

43 

84 

124 

3. 

Adam Krska 

G, Mikulov 4 

4 

6 

- 

1 

2 

11 

- 

28 

65 

86 

4. 

Benedikt Bares 

G Dobruska 8 

6 







14 

77 

69 

5. 

Vojtech Votruba 

G Jana Keplera, Praha 6 

- 

6 

- 

2 

- 

- 

1 

15 

56 

64 

6 . 

Vojtech Stransky 

G Brno, tr. Kpt. Jarose 2 

2 

6 

4 

- 

4 

10 

3 

31 

60 

62 

7. 

Adam flustava 

European School Luxern- 8 
bourg II 

— 

— 



2 

10 


20 

80 

56 

8 . 

Ales Opl 

Gymnazium Praha 3 6 

6 

6 






18 

87 

55 

9. 

Matej Charousek 

G Na Vitezne plani, Praha 2 

4 

2 

1 

1 

5 

4 

- 

19 

40 

43 


10 . Barbora Cemanovd 

11. Adela Kolembusova 

12 . Sara Byskova 

13. Filip Zikes 

14. Petr Sicho 

15. Eliska Durstova 

16. Milan Marek 

17. Jan Ptacek 

18. Zuzana Lisztwanova 
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Hranicna cena bananov, pri ktorej sa Mikulasovi tento podvod vyplati je 900CZK-kg _1 . Tato 
hranica je nezavisla 11 a mnozstve pouziteho lielia, cize ak by sine pouzili viac helia ako mame, 
vrecusko by sa zacalo vznasat. Zbytocnc by sme mrhali heliom a pani predavackc by sa to 
mohlo zdat podozrive. 


Bonus 

Z predchadajucich vypoctov vieme, ze musi plat it 


Cb 

C P 



> 1 , 


pricom Cp je cena plynu za kilogram a Qp je hustota plynu. Potrebujeme najst taky plyn, pre 
ktory sa podvod oplati pri cene bananov 30CZK*kg -1 . Najdeme si ceny 1 a hustoty 2 roznych 
technickych plynov a zistime, ci je niektory z nich vyhovujuci. Samozrejme, zaujimaju nas len 
plyny s nizsou hustotou ako vzduch q vm i = l,2759kg*m -3 , ktore su bezpecne na pouzitie v 
danej situacii. Ziskane udaje, z ktorych sme prepocitali cenu na jeden kilogram v CZK-kg 
a potom dosadili do vyssie uvedenych vztahov, umiestnime do prehladnej tabulky 1. 


Tab. 1: Tabulka plynov 


plyn 

£p 

Cp 

<?£ ( L>n,zd _ l\ 

Cp V up J 

kg-m -3 

CZK-kg" 1 

vodik 

0,08895 

6428 

0,06 

helium 

0,1762 

5100 

0,04 

zemny plyn 

0,7 

8 

3,2 

ace ty len 

1,147 

560 

0,006 

dusik 

1,234 

334 

0,003 


Plyn, pri ktorom sa to oplati, je - podia nami najdenych udajov o cenach plynov - zemny 
plyn, hlavne kvoli jeho nizkej cene. 


Katar'ina Castulikova 
katka.castulikovaQfykos.cz 


Ulolia III.2 . . . efektivni kafe 3 body; prumer 2,51; resilo 47 studentu 

Jsou dve hodiny v noci a Jachym si jde uvarit kafe. Na plotynku, kterou tvon litinovy valec 
o polomeru r a vysce h, polozi konvici s tcpclnou kapacitou Ck- Konvice obsahuje vodu o obje- 
mu V , ktera ma pocatccw teplotu T v . Zbytek soustavy ma pocatecni teplotu T s . Jakaje celkova 
ucinnost (tj. pomer encrgie pfijatc vodou ku dodane energii) ohfevu vody z jejf pocatecni tep- 
loty na teplotu T = 100 °C? Nezname hodnoty si dohlcdejte v tabulkach, ncbo je odhadnete. 
Pfedpokladejtc, ze dcj probcline tak ryclile, ze vsechny tcpclne ztraty muzeme zanedbat. Pro 
uplnost zadani necht T S ,T V < T Jachymovi dosel energy drink. 

1 Ceny plynu prevzaty z https://www.kurzy.cz/plyn/srovnani-cen a https://www.messer.sk/cennik 
2 llustoty prevzaty z http://www.convcrter.cz/tabulky/hustota-plynu.htm 
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Necht je hustota litiny #i, jeji merna tepelna kapacita je c\. Oelkova kapacita soustavy tak je 

C s = Ck 4- nr 1 2 ligici . 

Pfi ohrivani vody musime nutne ohrat i plotynku s konvici, k cemuz jim musime dodat teplo 
Q s = C S (T — T s ) = (C\ + nr 2 h 0i c\) (T - T s ) . 

Obdobne pro vodu mame g w a c w a musime ji dodat teplo 


Qv = V g v c v (T — T v ). 


Ucinnosti rozumime pomer energie, ktera se vyuzije na olifev vody, a celkove dodane energie, 
tcdy 


Qv 

Qv + Q s 


_P QyCy(r ly) _ 

Vg v c v (T - T v ) + (C k + nr 2 hg[Ci) (T - T s ) ' 


Nyni uz jen zbyva dosadit realisticke hodnoty. Nasc plotynka ma polomer r = 10 cm a vysku 
h = 2 cm. Konvice vazf 500g a je z oceli, jeji tepelna kapacita tedy je (7k = 230 J-K _1 . Na kafe 
potrebujeme V = 250 ml vody. Muzeme predpokladat, ze zacatku ma vse vcetne vody teplotu 
T v = T s = 20 °C. Mernou tepelnou kapacitu vody c v = 4 180.J-K -1 a litiny ci = 545J-K -1 , 
stejne tak jejicli hustoty g v = 1000kg-ni } a g\ = 7200kg-m~ 3 jsme nasli v tabulkach. 

Dosazenim hodnot do vzorce vyse dostavame odhad pro ucinnost ohrevu vody 77 = 0,28. 
Tato nadocekavani vysoka hodnota je zpusobena velkou tepelnou kapacitou vody v porovnani 
s kovy. 


Jachym Bdrtik 
tuakiQfykos.cz 


Ulolia III.3 ... toplicko v Dysonovc sforc 6 bodu; prumer 5,36; resilo 44 studentu 

Jaky polomer by rnusela nut Dysonov a sfera, a by obklopila hvezdu se zarivyin vykonem Slunce 
tak, ze na vnejshn povrchu teto sfery by byla teplota t = 25 °C? Neuvazujte pntomnost atmo- 
sfcry v Dysonovc sfcrc. Dysonova sfera by me la byt rclativnc tcnka duta struktura kulovcho 
tvaru obklopujici danou hvezdu. Karel ma rad Dysonovy sfery. 

Ze zadani plyne, ze teplota Dysonovy sfery se ustalila na konstantni teplote / = 25 °C. K vyre- 
seni ulohy nam tedy staci analyzovat tok energie Dysonovou sferou. Predpokladame, ze Dyso- 
nova sfera je schopna pohltit veskerou energii pfichazejici ze Slunce. Tato energie je prenasena 
v podobe elektromagnetickeho zareni. Ze zadani vime, ze dana hvezda ma zarivy vykon Slun¬ 
ce L = 3,83 • 10 2 ° W ? Tento vykon je rovnomerne vyzaren do vsech smeru. Mnozstvi energie, 
ktere dopadne za jednotku casu na jednotkovou j)lochu Dysonovy sfery je 

L _ L 
S 4 Tir 2 

kde S je plocha Dysonovy sfery a r je jeji polomer. Rozdil vnitrniho a vnejsiho poloineru 
zanedbavame. Dysonova sfera se navenek jevi jako absolutne cerne teleso, protoze veskerou 
energii z hvezdy pohlti a navenek vyzaruje pouze tepelne zareni. Podle Stefanova-Boltzmannova 

1 Viz napr. https://nssdc.gsfc.nasa.gov/planetary/factsheet/sunfact.html. 
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vach zostavit jedinu diferencialnu rovnicu prveho radu. Pre premennu r bude nasa vysledna 
diferencialna rovnica 

kde E je mechanicka energia a / moment hybnosti telesa - dva integraly pohybu. 

K rieseniu tejto rovnice sa dostaneme v serialovej ulohe a d’alsom diele serialu. Nasledne 
plynule prejdeme k problemu dvoch telies, rieseniu Keplerovej ulohy a mozno aj odvodeniu 
Keplerovych zakonov. 

Mnohyin z vas d’akujcm za reakcie 11 a vyzvu k spatncj vazbe, dojiny boli zatial celkom 
pozitivnc, za co som nid, pretoze je to pre mha povzbudzujuce. Samozrejme nad’alej plati, ze 
budem rad za vase navrhy, pripomienky a komentare k serialu. Som optimista a dufam, ze pocet 
riesitefov serialu sa nebude znizovat, tak si vsetci navzajom drzme palce. 
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je r je zrejme zjavne. Ostatne dva uhly si vieine jednoducho predstavit, ked si predstavime, 
ako to vyzera v planetariu (tain sa nad nami nachadza sfera). V tomto prfpade je polomer 
sfery pevny a nemeni sa. Ak stojimc v strode, tak spojnica nas a nejakeho vyznacneho bodu 
(v planetariu sa pouziva juh) je nasa kartezska os x. Suradnice nejakej hviezdy vieme urcit 
pomocou sferickych suradnic: hviezdu spojime najkratsou moznou ciarou po sfere s horizontom. 
Teda spravime jej priemet na „horizont“. Potom odcitame uhol medzi tymto jej priemetom 
a nami zavedenou osou x proti smeru hodinovych ruciciek, co bude nasa suradnica ip. Na zaver 
urcime suradnicu i) tak, ze od 90° odcitame uhol medzi horizontom a hviezdou po nami najdenej 
najkratsej spojnici. 

Snazil som sa to vysvetlit zrozumiterne a nazorne, na internete sa ale nachadza mnozstvo 
obrazkov pre tych z vas, ktori veci potrebuju vidiet, a nestaci im vizualizovat si ich. A mozno 
este jedna poznamocka k zavedeniu sferickych suradnic. Mozno vam prislo neintuitivne, ako 
som v poslednom kroku odcital uhol medzi hviezdou a obzorom od 90°. Samozrejme, dostaneme 
rovnako dobre pouzitelne suradnice, ak d definujeme priamo ako uhol medzi (pre nas pripad) 
horizontom a hviezdou. Prevodne vzfahy sa potom trochu pozmenia. Ak ste zvyknuti na tuto 
definiciu sferickych suradnic, pokojne ich pouzivajte, no pamatajte, ze nemusi byt na prvy 
pohl’ad jasne, ze vztaliy ziskane pri tychto roznych zavedeniach suradnic su totozne. 

Na zaver tejto casti si spomenieme este jeden priklad, z ktoreho vysledkov budeme intenzivne 
tazit nabuduce. 

Pohyb hmotnej castice v gravitacnom poli 

V roku 1889 Norsko-Svedsky krai Oskar II. vyhlasil k prilezitosti svojich 60-tych narodenin 
sutaz o najdenie analytickeho riesenia problemu troch telies vo forme konvergujuceho mocnin- 
neho radu, ktoru sam oznacil za (vo volnom jirckladc): ”Sutaz o dolezity objav v risi vysokcj 
matematickej analyzy”. Uz nasledujuci rok Henry Poincare publikoval clanok v casopise Ac¬ 
ta Mathematica, v ktorom dokazal, ze neexistuje dostatok integralov pohybu na analyticke 
vyriesenie problemu troch telies, co kral’a uspokojilo a odmena za ”vyhru” v sutazi mu bola 
udelena. 

Az v roku 1912 sa Finskemu matematikovi Karlovi Sundmanovi podarilo vyriesit tento 
problem, avsak formou radu s exponentmi, ktore su celociselnym nasobkom Tento rad vsak 
konverguje vel’mi pomaly, a aby dosahoval presnost porovnatelnu s numerickymi vypoctami, 
je nutne vziat priblizne jeho prvych iq 8()()0()()() clenov. Problem troch telies je tak prakticky 
riesitelny iba numericky. Existuju vsak zjednodusene pripady, ktorych riesenia s rozumnou 
presnostou mozeme ziskat aj analyticky. Takym je napriklad system Slnko-Zem-Mesiac. 

Nez sa ale dostaneme k tomuto problemu, musime si prejst cez problem dvoch telies. A pred- 
tym ako sa dostaneme k tomu, je potrebne si ukazat, ako vyzera Lagrangian hmotneho bodu, 
ktory sa pohybuje v nejakom potencialovorn poli. Budeme predpokladat, ze nase skumane teleso 
svojim pohybom toto potencialove pole nijako neovplyvhuje. 

Nase teleso sa bude pohybovat v centralnom poli s potencialom V (r) a zo skusenosti vieme, 
ze sa bude pohybovat v jednej rovine (tuto skusenost si v ulohach k serialu budete mat moznost 
overit). Na jeho popis preto pouzijeme polarne suradnice. V nicli bude jeho Lagrangian 

L = (r 2 + r 2 <p 2 ) -V(r). 

Tcraz uvidime silu integralov pohybu. Vo vseobecnosti by boli Lagrangeove rovnice v tomto 
pripade 2 rovnice druheho radu, kdezto ak urcime integraly pohybu, dokazeme z nich po upra- 
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zakona je celkova intenzita M tepelneho zareni cerneho telesa, coz je celkova energie vyzarena 
za jednotku casu jednotkovou plochou zdroje zareni, rovna 

M = oT 4 , 


kde <t = 5,07 • 10 -8 W-m -2 -K -4 je Stefanova-Boltzinannova konstanta a T = 298 K je termo- 
dynamicka teplota Dysonovy sfery. 

Zarivy vykon hvezdy a vykon tepelneho zareni Dysonovy sfery jsou jedine zdroje zareni, 
ktere musime zapocitat. Prestoze princip fungovani Dysonovy sfery je schovan v poinyslne cerne 
skrihee, staci nam vedet, ze Dysonova sfera je schopna pohltit veskerou energii ze Slunce. Na 
druhou stranu, aby byl nas vypocet korektni, musi tepelne zareni vnejsi plochy Dysonovy sfery 
navzdy opustit tuto soustavu, coz jsme mlcky predpokladali. Nas prcdpoklad byl v poradku, 
nebot koule je konvexni teleso a energie vyzarovana malou ploskou AS na jejim povrehu, sfere, 
unika do prazdneho poloprostoru. To nam zarucuje, ze Dysonova sfera neozafuje sama sebe 
vnejsim povrehem. 

Nyni konecne dejme do rovnosti energii z hvezdy, ktera dopadne za jednotku casu na jed¬ 
notkovou plochu Dysonovy sfery, a intenzitu M tepelneho zareni Dysonovy sfery, 



Odtud ziskame hledany polomer r Dysonovy sfery, 



Po ciselnem dosazeni dostavame, ze hledany polomer Dysonovy sfery pro zadane podminky 
je r = 2,61 • 10 11 m = 1,74 AU. 

Vaclav Mikeska 
v.mikeskaOfykos.cz 


Ulolia III.4 ... destrukce smycky 8 bodu; prumer 6,00; resilo 25 studentu 

Predstavme si medenou smycku o polomeru r, ktera jc urcen a rovinou, na mz je kolme inagne- 
ticke pole s magnetickou indukef B. Maximalm povolene tahove napeti ve smycce je a p . Nyin 
zacneme menit magncticky tok vc smycce z puvodm hodnoty <I>o podle vzalm 4>(^) = 4>o + at, 
kde a jc klndna konstanta. Urcete, za jak dloulio dosahneme vc smycce niaxirnalm'ho tahovcho 
napeti. 

Napoveda: Napeiovou si'lu vc smycce muzeme spocitat jako T = \BIr\. 

Vitek vzpomind na A P Physics. 

Zamysleme se nad tim, co se v uloze deje. V okamziku, kdy zacneme menit magneticky tok 
prochazejici smyckou, se v ni zacne indukovat napeti a tedy i proud. Tuto skutecnost nam 
popisuje Lenzuv zakon, ktery rika, ze indukovany elektricky proud v uzavrenem obvodu ina 
takovy smer, ze svym magnetickym polem pusobi proti zmene magnetickeho indukeniho toku, 
ktera je jeho pricinou. System se snazi jakoby vratit do puvodniho stavu, a proto rna Lenzovska 
indukce smer opaeny iiuvodnimu B. Reel matematiky 


£ 


d $(0 

d t 
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Indukovany proud ve smycce pak urcfme snadno pomocf Ohmova zakona. Oznacfme-li odpor 
smycky jako /?, pak 

j _ Jfl _ <*A 

R Q‘2nr ’ 

kde A je prurez smycky a g rezistivita materialu. Je zrejme, ze T zavisf na B , a to zavisf na 
case. Pro moment to, kdy dosahneme maximalmho tahoveho napetf, musf platit nasledujfcf 
podmfnka 

T(to) B{to)I r B(t 0 )a 

Gp A A 2 qk * 

Dale musfme urcit B(to ), z definice 4 zfskame 


B(to) = 


*{to) _ *0 + cxto 


ti r* nr* 

Nynf stacf dosadit za B(to) do vztahu pro tahove napetf, to pak prejde na 


(4>o + ato)a 
Gp ~ 2qk 2 r 2 ’ 

Cflem uloliy bylo zjistit, za jaky cas se toto stane. Proto ze vztahu vyse vyjadrfme to, finalnf 
vyraz pak je 

2t i 2 gr 2 a p <I>o 

to = -o-• 

or a 

V rcscnf jsmc mlcky zanedbali vlastnf magneticke pole indukovaneho proudu. Z prurezu A 
bychom si mohli dopocftat vlastnf polomer vodice a odhadnout velikost intenzity vlastnfho 
magneticke pole, nicmene typicka smycka ma vlastnf polomer mnohem mensf nez polomer 
R a tudfz muzeme tuto skutecnost zanedbat. Pri resenf jsrne taktez zanedbali vlastnf zmenu 
rozmeru smycky. 


Pro zajemce 

Tyto radky jsou venovane resitelum, kterf by radi vedeli, jak si odvodit napovedu v zadanf. 
Jak jiz vfme z textu vyse, v uzavrene smycce se zacne indukovat proud. Proto zacne na vodic 
pusobit sfla ve srneru do stredu smycky (viz. obrazek) podle Amperova zakona sfla dana jako 

dF = I d/ x B, 

kde d/ je maly element smycky. Ten si muzeme vyjadrit pomocf malcho useku d0 jako rd0. 
Protoze jsou na sebe vektory v rovnici kolme, tak se vztah zjednodusf na 

dF = I BrdG. 

Z obrazku nynf muzeme vyjadrit tahovou sflu ve smycce v zavislosti na malem uhlu d0. Hori- 
zontalnf slozky se vyrusf (jsou rovny T cosd0/2) a vysledna sfla dovnitr smycky (pusobfef proti 
dF) je dana souctem dvou vertikalnfch slozek, zfskame tedy 

dr = 2T sin ^ « Td0 , 

4 Magneticky indukeni tok vytvafcny magnctickou indukef B na libovolnc orientovane plose S jc definovan 

jako $ = B ■ dS. Pro homogenni pole a rovinou plochu plati 4> = BS cos a (uhel or svira normalovy vektor 
J s 

plochy s vektorem magneticke indukce). V nasem pripade jsou na sebe vektory kolme. 
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Vidfme, ze tento vyraz sa bude stale zachovavat, nakol’ko je po celej trajektorii konstantny. Na 
zaver mozno este poznamka k nazvu „zovseobecnena hybnost“. Kedzc Lagrangian ma fyzikalny 
rozmer energie (jedna sa o rozdiel dvoch energif), ked’ zderivujeme „energiu“ podia „rychlosti“, 
dostaneme nieco, co ma fyzikalny rozmer hybnosti. 

Okrem toho existuje este jeden integral pohybu? 0 a ten sa nazyva „zovseobecnena energia“. 
Zaujfmavejsie nez jeho konstrukeia je, ze k Lagrangianu existuje zachovavajuca sa zovseobec- 
nena energia prave vtedy, ked Lagrangian explicitne nezavisf na case (teda na case zavisi len 
prostrednfctvom zovseobecnenych suradnfc a rychlostf). Jej definfeia je 

h{q l ,q i ) = Y j ^~q l -L. 

1=1 

Znova si vsimnene, ze tato velicina ma fyzikalny rozmer energie, takze jej nazov dava dobry 
zmysel. 

Este predtym, ako realne zacneme riesit prfklady, je podia mna dolezite ukazaf (pre niek- 
torych zopakovat) novy sposob hladania zovseobecnenych suradnfc. 


Najcastejsie krivociare su rad nice 

V mnohych ulohach sa stretavame s pohybom po povrehu sfery alebo po kruznici, alcbo s pro- 
blemom, ktory ma rotaenu symetriu. Hlboko v podstate prfrody je totizto zakorenena jedna 
vlastnost, ktorej sa snad’ budeme venovat v poslednom diele serialu a z ktorej sa da odvodit 
cela fyzika. Zjednodusene by sa dalo povedat, ze prfroda je „leniva“, a preto necha planety obie- 
hat okolo Slnka po elipsach a nie po naprfklad stvorcoch. Vd’aka tejto podivuhodnej vlastnosti 
sa zakriveny pohyb najcastejsie odohrava po elipsach velmi blfzkych kruzniciam (mimo me- 
chaniky na Zemskom povrehu, kde uvazujeme liomogenne gravitacne pole). Preto je castokrat 
vyhodnejsie poeftat v suradniciach, ktore dobre vystihuju symetriu danej sfery alebo kruznice. 
Zavadzaju sa preto polarne a sfericke suradnice. 

Polarne suradnice definujeme 


x = r cos(<p), 
y = r sin(<p). 

Ak si v rovine nakreslfme kartezsku mriezku, potom bod s kartezskyrni suradnicami [x,y\ 
mozeme popfsat aj pomocou uhlu ip zovreteho mcdzi kladnym smerom osy x a polpriamkou 
z pociatku do nasho bodu; uhol mcriame proti smeru hodinovych rucicick (smerom od osi x 
k osi y). Suradnica r je potom vzdialenosf tohto bodu od nuly merana po ich spojnici. 

Druhe, o cosi tazsie na predstavu, su sfericke suradnice. Prevod mcdzi sferickymi a kartez- 
skymi suradnicami je nasledovny 


x = r cos ((p) sin('^), 
y = rsin(cp) sin(^), 
z = rcos^-d) . 

Mozeme si vsimnut, ze je to podobne ako pri polarnych suradniciach, len este prenasobene 
sfnusom a kosmusom d’alsieho uhla, ktory sine si oznacili i). To, ze vzdialenost bodu od pociatku 

30 V skutocnosti mozu existovat aj ine integraly pohybu, spomenute sa vsak daju najst pomerne jednoducho. 
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Rychlost dopocitame tak, ze prederivujeme rovnicu pre polohu kyvadla, aby sine dostali 
rovnicu pro uhlovu rychlost kyvadla 

= c2 \/? cos (/?') - vT v ° sin (/fO • 

Z tejto rovnice a podmienky, ze rychlost ma byt na pociatku nulova, dopocitame druhu kon- 
stantu. Cely druhy clen rovnice bude po dosadeni t = 0 rovny nule, zatialco prvy bude nulovy 
len ak C 2 = 0. Nase riesenie pre dane pociatocne podmienky bude teda vyzerat 

<p(t) = fOCOS (y|t) ■ 

Prave sme si podrobne ukazali jeden zo sposobov riesenia Lagrangeovych rovnic. Jedna 
sa o to najjednoduchsie, co nam rnoze napadnut. Hlavnou nevyhodou je ale to, ze riesenie, 
ktore takto najdeme, nie je presnym riesenim, a spravne vysledky dava len pre body blizke nule 
(pretoze spravidla linearizujeme funkciu v nule). Existuje ale d’alsi trik pouzitelny pri niektorych 
prikladoch, ktory nam vyraznc ulahcuje ich analyticke riesenie. 


Integraly pohybu 

Jeden z trikov, ktore mozeme pri rieseni Lagrangeovych rovnic pouzit, je ziskat dodatocne 
rovnice, ktore nam zjednodusia riesenie samotnych Lagrangeovych rovnic. Zo strednej skoly 
viete, ze pri rieseni uloh z mechaniky ste castokrat vyuzivali zakony zachovania. Univerzalny sa 
moze zdat zakon zachovania mechanickej energie, ktory ale neplati vzdy. Dalej su uzitocne aj 
zakony zachovania hybnosti a momentu hybnosti. Castokrat nebolo jednoduche zistit, najma 
pri vehni komplexnych problemoch, ktora zlozka hybnosti sa zachovava a preco. Pripadne ci sa 
zachovava alebo nezachovava energia. V Lagrangeovom formalizme je tento problem nalezite 
osetreny a z matematickych vlastnosti Lagrangianu (Lagrangeovcj funkcie) vieirie jednoducho 
urcit, ktore veliciny sa zachovavaju pocas celeho dcja. Tieto veliciny nazyvame integraly pohy¬ 
bu. Matematickejsie povedane, integral pohybu hladame ako funkciu }{q 3 ,4j) zovseobecnenych 
suradnic a zovseobecnenych rychlosti, pricom od tejto funkcie pozadujeme, ze na skutocnej tra- 
jektorii riesiacej Lagrangeove rovnice bude nadobudat po cely cas jednu konstantnu hodnotu. 
Povedali sme si, ze v praxi sa jedna napriklad o hybnost telesa. Z Lagrangeovych rovnic vidime, 
kedy sa nieco ako „zovseobecnena hybnost“ zachovava, celkom rychlo. 

Predpokladajme, ze Lagrangeova funkcia zavisi na vsetkych zovseobecnenych suradniciach 
okrem nejakej jednej q r . Potom z-tu Lagrangeovu rovnicu ziskame 


d_^L _ dL_ 
d t d(jr dqi 


= 0. 


Derivacia Lagrangianu podia tejto suradnice bude teda 0 (pretoze pri zmene tejto suradnice sa 
Lagrangian nezmeni). Z toho dostavame 


T^=»- 

d t dqi 


Na zaver preintegrujeme podia casu 


5L 

dqi 


= konst. 
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kde jsme pouzili aproximaci pro male hodnoty argumentu funkce sinus. Z podmmky rovnovahy 
j)ak mame 

I BrdQ = TdQ => T — I Br . 

A mame dokazano. Tento princip je ve fyzice velmi casty a doporucuji ho radne prostudovat. 0 


T 



Vit Beran 
vit.beran@fykos.cz 

Uloha III.5 . .. bodova 8 bodu; prumer 5,32; resilo 22 studentu 

Uvazujmc hmotny bod umisteny v jednodimenzionalnim prostoru. Jcho pocatecru pozice i rycli- 
lost je nulova. Bod se dokaze pohybovat s libovolnym zrychlemm z intervalu [—a , a]. Nazveme 
M{t) mnozinu vscch moznych stavu (x,v) takovych, ze bod se v case t muze nachazet na pozi- 
ci x s rychlosti v. Sestrojme graf zavislosti v na x v case t. Mnozina M(t) v tornto grafu vytvorf 
plochu S(t). Analyticky popiste krivky ohranicujici S(t). 

Bonus Najdctc funkciu zavislost obsahu S(t). 

Jdchym chtel jistou trividlm ulohu resit jako specidlm pfipad teto. 

Nejdrive podotkneme, ze pokud je mozne nejakeho bodu (x,u) dosahnout v case i, zrejme 
je mozne ho dosahnout i v jakemkoli vetsim case tim, ze se pred zacatkem pohybu budeme 
odpovidajici casovy rozdil pohybovat s nulovym zrychlenim. Z toho vyplyva, ze mnozina M se 
v case pouze zvetsuje. 

Zvolme nejaky konkretni cas t. Nyni si predstavme, ze se po celou dobu od pocatku do 
casu t pohybujeme pouze se zrychlenim a. Neni tezke spocitat, ze dosahneme rychlosti u m ax = 
= at a polohy 

— ^ f 2 

^max — fy dt . 

Toto je zrejme nejvetsi rychlost a zaroven nejvetsi vzdalenost, ktere muzeme v danem case 
dosahnout, takze musi byt soucasti okraje mnoziny. Obracene, nejdale v zapornem smeru se 
muzeme dostat do bodu ve vzdalenosti x m i n = —x ma . x s rychlosti v m in = —v ma . x . 

Dale zkusrne najit, jakou nejvetsi rychlost muzeme mit na souradnici x. Pro kazdy bod 
od x m in do x max tak najdeme liorni hranici mnoziny pro rychlost bodu. Situace je zrejme 

5 Dalsim typickym prikladem na nej je napr. capstan equation. 
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symetricka vuci soucasne zmene znamenka u rychlosti a polohy, takze dolni hranici rychlosti 
najdeme pouhym bodovym zrcadlenim vysledku vzhledem k pocatku. 

Abychom pfi pohybu z klidu dosahli ncjake rychlosti, musime zrychlovat a pri tom nutne 
urazime nejakou vzdalenost. Muzeme si rozmyslet, ze pokud bychorn se pohybovali s mensim nez 
maximalnim zrychlenim, jakekoli rychlosti bychom dosahli na vets! vzdalenosti nez v pripade 
pohybu s maximalnim zrychlenim. Jinak receno, chceme-li na co nejkratsi vzdalenosti dosahnout 
co nejvetsi rychlosti, musime se pohybovat s maximalnim zrychlenim. 

Snadno spocitame, ze pro zrychleni z nuly na v musime urazit vzdalenost 


Jakou maximalm rychlost tak muzemc mit napriklad v pocatku, tedy v bode x = 0? Jednoduse 
musime nejdrive zrychlovat na opacnou stranu, cimz se dostaneme do zapornych hodnot osy x. 
Potom muzeme zacit zrychlovat v kladnem smeru, cimz se postupne vratime zpet do pocatku, 
tentokrat uz s kladnou rychlosti v. 

Nyni zase vyvstava otazka, jak se v co nejkratsim case dostat co nejdale od pocatku, 
abychom se pak mohli vratit s co nejvetsi rychlosti? Odpovecf je stejne prosta, jako cely zby- 
tek ulohy - pohybem nejdrive s maximalnim zapornym zrychlenim a lined pote s maximalnim 
kladnym zrychlenim. 

Pojd’me konecne uvahy vyse prcpsat do rovnic. Pro daily cas t si zvolime bod x a maximalm 
rychlost, kterou v nem muzemc mit, oznacime v. Po cas t\ budeme mit zrychleni —a a po dalsi 
cas £1 zrychleni a. Tim dosahneme bodu (y, 0). Odtud se budeme pohybovat se zrychlenim a 
po dobu t 2 a konecne dorazime do bodu (x,v). Dostavame rovnice 

y — ~2 a ^ + ~ a + 2 ^ 1 ) = ~ a ^i » 

x-y= ^at 2 , 

v = at 2 , 
t = 2t!+t 2 . 


Znamc a, i a .t, nezname u, £ 1 , t 2 a y. Ctyri rovnice pro ctyri neznamc jsou presne tolik co 
potrebujeme, abychom si z nich vyjadrili v. Nicmene zjistime, ze se neobejdeme bez odmocnin. 
Nam ale staci najit vztah pro hranicni body v case £, teda zbavit sa pomocnych neznamych £ 1 , 
t 2 a y. Dostat muzeme napriklad rovnici 

x = - 7 - (v 2 + 2 vat — a 2 t 2 ) . 

4a v 7 

Coz je samozrejme parabola. Jeji osa je rovnobezna s osou x, kterou parabola protina v bo¬ 
de (xi, 0 ), kde 

1 ,2 

Xl = --at , 

coz je nejvzdalenejsi bod v zapornem smeru, kam se muzeme dostat s nulovou konecnou rych¬ 
losti. Odtud vede hranice az do bodu (.T max , ^max), ktery jsme urcili drive. Parabola samozrejme 
pokracuje i dal, tarn vsak nerna fyzikalni smysl (prislusny cas 1 1 by vysel zaporny). Nicmene 
stale jeste jsme nenalezli druhou hranici mnoziny od bodu (:r max , Vmax) zpet na osu x. 
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Matematicke kyvadlo - In medias res 

V minulej casti sme si ukazovali, ako postupne najst zovseobecnene suradnice, z nich cisto 
mechanicky zostavit Lagrangeovu funkciu a z nej postupnym derivovanim zostavit Lagrangeove 
rovnice. To sine si ale neukazali na nasom priklade matematickeho kyvadla. Dokoncime to teraz. 
Pre pripomcnutie, Lagrangeova funkcia kyvadla mala tvar 


1 2 2 

L = —ml ip + rrujl cos (ip) 
a Lagrangeove rovnice zostavujeme nasledovne 


d_dL _ DL _ 
c \t dq d(j 

Dosadime za L Lagrangeovu funkciu matematickeho kyvadla a po prederivovani dostaneme 

ml 2 (p + myl sin(c^) = 0 . 


Vidiine, ze sa jedna o obycajmi diferencialnu rovnicu druhcho radu, prctoze najvyssia derivacia 
nasej hladanej funkcie p(t) je druha derivacia. Co je ale ovel’a vacsi problem, rovnica je zaro- 
veii aj nelinearna , a to prave jej druliy clen, ktory je tvoreny sinusom nami hladanej funkcie. 
Standardne sa nedaju podobne diferencialne rovnice riesitl analyticky, t.j. pomocou znamych 
dobre definovanych funkcii. R.iesit ich mozme napriklad numericky s presnostou, aku uzmime 
za vhodnu. Inou moznostou je riesit rovnicu analyticky tak, ze pred tym, ako sa pustime do 
riesenia, danu nelineanu 29 funkciu aproximujeme. Toto je okamih, ked’ prvykrat v praxi po- 
uzijeme Taylorov rozvoj. Rozvinieine nasu funkciu do Taylorovho radu. Vieme, ze cim sa od 
bodu, v ktorom funkciu rozvijame, nachadzame d’alej, tym viac clenov jej rozvoja potrebujeme 
pre dosiahnutie pozadovanej presnosti. Naopak, ked sme blizko stredu tolito rozvoja, staci nam 
clenov mcnej. V istej malej vzdialenosti nam dava dostatocnu presnost (cim vyssiu presnosti 
pozadujemc, tym monsie je okolie bodu, kde je aproximacia presna) aj prvy clen taylorovho 
rozvoja. Ten je obvykle linearny, pripadne konstantny (v pripade konstantneho clena vezmcme 
prve dva cleny). V tomto pripade vezmeme z rozvoja sinusu len prvy linearny clen. Dostavame 
linearizovanu diferencialnu rovnicu 


ml 2 (p + my Ip = 0 . 

Tato rovnica sa da jednoducho riesit. Na riesenie diferencialnych rovnic sa vsak v tomto seriali 
nebudeme sustredit, ale nechame ich riesenie na pocitac. Osobne preferujem, ako ste si uz mohli 
vsimnut, na riesenie rovnic Wolfram Alpha, ktory nam v tomto pripade da vysledok 

<p(l) = c 2 sin cos • 

Kde c*i a C 2 su integracne konstanty. Ich hodnoty si vieme jednoducho dopocitaf. Napriklad tak, 
ze budeme pozadovat, aby na zaciatku rnalo kyvalo vychylku po a uhlovu (a teda aj klasicku) 
rychlost p = 0. Z prvej podmienky vidime, ze ci = c^o, co dostaneme po dosadeni t = 0 do 
nasej pohybovej rovnice. 

29 Medzi nclinearne funkcie samozrejme patria aj kvadraticke, kubicke a ine mocninne funkcie. Nemoznosf 
analyticky riesit rovnice je vsak najma pri nclincarnych rovniciach, kde nclincarita je sposobena inou ako 
polynomialnou funkeiou. 


31 





Fyzikalni korespondencm seminar MFF UK 


rocnik XXXII cislo 4/7 


polynom pozostavajuci napriklad z piatich clenov 


x 3 x 5 x 7 X 9 
sm(x) ~x-- + — - - + 


Nie je tazke si vsimnut vzor v spravani tohto polynomu. Ak by sme cliceli aproximovat tuto 
funkciu na celorn intervale presne, mohli by sme pokracovat v pridavani clenov polynomu az 
do nekonecna. Ukazuje sa, ze takyto nekonecny polynom sa naozaj bude rovnat funkcii sinus 
dokonca na celom jej definicnom obore, teda vsetkych realnych cislach. V ’’modernej mate- 
matike” je prave takyto nekonecny polynom (rad, sucet postupnosti) definiciou funkcie sinus 
(spolu s inymi definiciami, ktore su tejto ekvivalentne). Sinus vieme potom zapisat pomocou 
nekonecneho Taylorovho radu nasledovne 


sin(x) = ^(-l) n 

n=0 


x 2n+1 
(2 n+ 1 )! ' 


Pre usilovneho citatela zadavam ulohu rozmyslief, si, najlepsie rozpisat na papier, ze prvych 
5 clenov radu je naozaj rovnakych ako su cleny v polynome aproximujucom sinus uvedenom 
skor. Dalej mozem taktiez odporucit usilovnemu citatelbvi, aby si vycislil liore uvedeny polynom 
v bode n a zaroven si vycislil v bode 7 t aj sinus. Budete mat lepsiu predstavu, ako dobre Taylorov 
rad aproximuje sinus. 

Usilovni citatelia mi urcite mozu potvrdit, ze vysledok bol v oboch pripadoch vefmi podobny. 
Pre sinus to bola nula, pre polynom to neprezradim, ale prezradim, ze koren tohto polynomu 
je ~ 3,1487, co sa k cislu 7 t blizi vehni posobivo. 

Castokrat nie je jednoduche zistit Taylorov rad trikom alebo pozorovanim (aj v tych malo 
pripadoch, ked’ sa to da, to dice velmi verku skusenost a vela napocitanych prikladov). Existuje 
ale mctoda, ako vieme postupne dopocitavat jeden za druhym cleny Taylovovho radu. Uvedieme 
rovno jeho vzorec a jednym dychorn dodavame, ze si mozete vsimnut, ze dany rad je zvoleny 
tak, aby mal vsetky derivacie v danom bode rovnake ako povodna funkcia. 

Zobcrme Fubovolnu funkciu f(x) v bode (x = a), ktora ma v tomto bode vsetky derivacie 
a vsetky su konecne (dovod je jasny zo samotneho vztahu - ak by derivacie neexistovali, nemal 
by zmysel). Na intervale obsahujucom bod a je potom Taylorov rad tejto funkcie 


7/(x) = ^/< n) (a) 

n=0 


{x - a) n 
n! 


kde f^ n \(i) je n-ta derivacia nasej funkcie v bode a. 

Kazdu normalnu (tzv. analyticku) funkciu / vieme potom zapisat ako /(:/;) = Tf(x). Exis- 
tuju aj funkcie, pre ktore to neplati, ale s nimi sa nestretneme. Overte si, prosim, ze ak tento 
vzorec aplikujeme na sinus v bode a = 0 , dostancme postupne polynom, ktory sme uviedli ako 
aproximaciu sinu. Pouzime tento vzorec este pre odvodenie Taylorovho radu pre exponencialu 
e x v bode a = 0. Derivacia eponencialy je exponcnciala. Exponcnciala v bode 0 je rovna jcdnej. 
Taylorov rad exponencialy bude teda vyzeraf nasledovne 


X 



71=0 


Ako cvicenie odporucam odvodit si zo vzorceka taylorov rad pre cos(x) a In (1 + x) so stredom 
v bode a = 0 . 
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Pokud se budeme pohybovat polovinu casu t se zrychlenim a a druhou polovinu se zrychle- 
niin —a, skoncimc opet na ose x, tedy s nulovou rychlosti, na souradnici 





Bod (x2,0) tak zrejmc prcdstavujc ncjvctsi vzdalcnost, kam sc muzcmc od pocatku na ose x 
dostat. Ilranice mnoziny musi vest z bodu (x ma x,^max) do bodu (x 2 , 0 ). 

Nyni stojime prcd podobnym problemem jako v prvni casti reseni - snazime se najit nejmensi 
rychlost, kterou muzeme mit pro dane x mezi X 2 a .T max . Jednodussi vsak bude opacny postup 
- zkusme pro danou rychlost mezi 0 a u max najit nejvetsi x, do ktereho se muzeme v case i 
dostat. 

Jiste je pravda, ze maximalniho x dosahneme tehdy, pokud se v kazdem bode trajektorie 
budeme pohybovat s nejvetsi rychlosti, s jakou to bude prave mozne. Zacneme tedy zrychlovat 
se zrychlenim a, a to az do casu r. V case r naopak zacneme brzdit se zrychlenim —a, takze 
v case t budeme mit rychlost v. Prvni cast drahy se pohybujeme s maximalnim zrychlenim 
a tedy i maximalm moznou rychlosti. Naopak ve druhe casti brzdime nejvice jak to jde, a tedy 
kdybychom meli jeste o neco vetsi rychlost, uz bychoin to do casu t nedokazali ubrzdit na 
rychlost v. 

.lednoduchou uvahou jsmc ukazali, ze tento zpusob pohybu nam pro danou koncovou rych¬ 
lost v zajisti nejvetsi moznou urazenou vzdalcnost. Opet dostavame sadu rovnic 

x = iar 2 + (^ar (t - r) - (t - r)“ j , 
v = ar — a (t — r) = a (2r — t) . 


Nyni mame jen dve nezname, a sice x a r. Muzeme si tak vyjadrit 

x = — ( a 2 t 2 + 2 vat — u 2 ) . 

4 a x ' 

Pokud nam tento vysledek prijde povedomy, tak zcela opravnene. Jedna se o parabolu, ktera 
nam vysla jiz v prvni casti, bodove prevracenou podle pocatku. Jak jsme zminili vyse, ze 
symetrie situace vyplyva, ze tato parabola tvori dolni hranici mnoziny od bodu x m in do bodu 
X '2 , zatimco prvni parabola vytvari dolni hranici od bodu x m i n do bodu x\. 

Abychom shrnuli predchozi uvahy, mnozina vsech bodu, ve kterych se muzeme v case t na- 
cliazet, je ohranicena dvema parabolami, bodove symetrickymi vuci pocatku. Jejich analyticky 
predpis je 

x = y- ( v‘ 2 -f 2 vat — a 2 t 2 ) , 

4a v ' 

x = — (a 2 t 2 + 2 vat — u 2 ) . 

4a v 

Pro uplnost dodejme, ze pro kazdy bod uvnitr ohraniceneho prostoru muzeme najit takovy 
cas r, ze r < t a zaroven dany bod lezi na jedne z hranicnich parabol pro cas r. r Io je podle 
vubec prvni poznamky tohoto reseni dukazem, ze cela plocha nema zadne vnitrni hranice kolem 
oblasti, ktere by do hledane mnoziny nepatrily, ale zaroven by byly uvnitr vnejsich hranic. 
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Obr. 2: Ilranice vysledne mnoziny M v nejakem case t. 

Bonus 

Pokud jste to docetli az sem, tak jiste tusite, ze bonus je jiz jen trivialni zalezitosti. Staci vhodne 
zvolit podle ceho integrovat, abychom si zbytecne nepridelavali praci. Definujme prvni parabolu 
jako x = fi(v) a druhou jako x = h{ v )' Potom pro obsah plochy S zfejme platf 

r v m&x -| r v max i 

S = 2 {f 2 (v) - fi{v))dv = - (a 2 t 2 - v 2 ) du = - 

Jo a Jo a 

Jachyrn Bdrtik 
tuaki@fykos.cz 


2,2 v 
a t v —— 


2 2,3 

= —a t . 
o 


Ulolia III.P . .. osobin powcrbanka 10 bodii; prumcr 4,04; resilo 27 studentu 

Posledm procenta baterky v mobilu dochazi, powerbanku mate vybitou nebojste si ji pro jistotu 
nechali doina a 230 take neni nikdc v dohledu. Nebylo by skvclc nut ncustale pfi sobc vlastm 
zdroj clcktricke cncrgic? 

• Navrhnete nekolik ruznych zanzeni, ktera by dokazala vyrabet elektrickou encrgii pouze 
ze zdroju vascho tela. 

• Diskutujte jejich maximalru vykon a ucinnost. Co vsechno byste s jejich pomoci dokazali 
zasobovat elektrinou? 

• Diskutujte jejich dopad na vase zdrav! a fyzickou kondici. Kterc organy by vam v duslcdku 
jejich pfetezovam selhaly nejdnve? 
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Serial: Mechanika 


V tejto casti serialu dokoncime priklad, ktory sme minule zacali — vypocet matematickeho 
kyvadla. K tomu ale budeme potrebovat vediet, co je to Taylorov rozvoj. Dalej si ukazeme, 
ako sa daju s pomocou Taylorovho rozvoja vyriesit ulohy, ktore analyticky riesitelne nie su. 
Potom sa pozrieme na to, ake nenumcrickc mctody riescnia Lagrangcovych rovnic existuju, 
a zacncme sa venovat pouzitiu teoretickej mechaniky v nebcskcj mechanike, co bola historicky 
hlavna motivacia, kvoli ktorej bol Lagrangeov formalizmus vyvijany. 

Taylorov rozvoj 

Toto bude trochu matematicka vsuvka do serialu, kde sa citatel’om, co vedia derivovat, pokusim 
intuitivne vysvetlit, co to je Taylorov rozvoj. Tiez ukazem, ako sa pocita, ale najma uvediem 
zakladne vztaliy pre Taylorov rozvoj, ktore sa v teoretickej mechanike zidu pri pocitani prikla- 
dov. 

Asi najintuitivnejsi sposob, ako chapat a predstavovat si funkcie, je predstavit si ich graf. 
Z grafu clovck jasne vidf, ako sa funkcia sprava, kde klesa. alebo stupa, ako vel’mi klesa alebo 
stupa, ake ma korcne a podobne. Mnohe funkcie su ale velmi narocne na analyzu a vypocet ich 
presnych hodnot, pripadne sa velmi narocne integruju (co nam znemozhuje, ci vyrazne stazuje 
riesenie diferencialnych rovnic, ktore obsahuju taketo funkcie). Na druhej strane, existuje trieda 
funkcii, ktore sa daju velmi trivialne integrovat a derivovat - polynomy (mnohocleny). S tymito 
funkciami ste sa uz vsetci stretli na konci zakladnej alebo zaciatku strednej skoly. Vlastnosti 
a rozne manipulacie s linearnymi a kvadratickymi funkciami su vam teda urcite vehni dobre 
zname a netreba vas presviedcat, ze vypocitat hodnotu rubovofneho polynomu v danom bode 
bez kalkulacky vie s fubovolnou presnostou kazdy z viLs. Keby som vam dal ale za ulohu zratat 
na papieri sin(0,82), nikto (kto nepozna Taylorov rozvoj) by to nijako nedokazal. 

Mozno ste sa niekedy hrali s nejakym programom na vykresl’ovanie grafov funkcii a vsimli 
si, ze mnohe funkcie maju podobny tvar. Exponenciala je napriklad funkcia, ktora velmi rychlo 
rastie. Aj x 2 , x A alebo ine vyssie stupne polynomov rychlo rastu. Ak si napriklad cxponencialu 
v takomto programe (napriklad WolframAlpha) vykreslite a vo vcdlajsom okne sa budete hraf so 
scitavanim polynomov roznych stupnov, s trochou sikovnosti a invcstovaneho casu sa dostanete 
k polynomu, ktoreho graf sa bude velmi podobat na graf exponencialy (aspon vo vykreslenej 
oblasti, pre vefke x exponenciala zacne rast ovda rychlejsie ako akykol’vek polynom). 

Tento fakt si l’udia vsimli uz ovel’a skor. Na takejto podobnosti spravne zvolenych poly¬ 
nomov a dancj funkcie je postaveny Taylorov rozvoj. V skratke sa jedna o to, ze rubovol’nu 
(dostatocne slusnu) funkciu vieme na nejakom intervale lubovofne dobre aproximovat urcitym 
polynomom. Tato vlasnost je veTmi uzitocna, pretoze castokrat nas funkcie zaujimaju len na 
nejakych kratkych intervaloch. Napriklad sinus. Sinus ma periodu, takze pre 1’ubovoFne vel’ke 
cislo ho vieme spocitat, ak dokazemc spocitat sinus vsetkych cisel na intervale o dlzke periody 
sinu 2 k (napr. (—tc, k)). Sinus by nam teda teoreticky stacilo dobre ai>roximovat len na tomto 
kratkom intervale. V pripade sinusu, ktory chceme aproximovat v okoli nuly, rnozeme pouzit 
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Jciko jedno z moznych zanzeni uvazujte soustavu drobnych turbm umfstenych v kreviu'ni 
recisti. Vscchny argumcnty podpoftc co ncjpfesnejsimi vypocty. 

Jachym mcl pocAt, ze mu chybi nejakd energie. 

Lidske tclo ziskava energii z preineny chemickych latek v potravinach na jine. Tuto energii opet 
uchovava ve forme chemickych vazeb. Ncjvetsi cast energie (kolem 80 %) se nakoncc premcni na 
teplo a zbytek se spotrebuje pri vyparovani vody? Stalo by za to toto teplo vyuzit a premenit 
ho na elektrickou energii. Dale se nabizi pouzit vhodne umistene turbinky (napriklad v krevnim 
recisti nebo v dychaci trubici). Energii muzeme ziskat i spalenim methanu a vodiku, produko- 
vanych ve strevech. Jednu metodu jiz dlouhou dobu vyuzivame a tou je mechanicka energie 
neboli prace svalu - urcite jste si jiz nekdy tocili klikou nabijejici svitilnu. Nektere z techto 
napadu si detailneji rozebereme. 

Nejdfive se podivejme na nektera zarizeni, ktera bychoin chteli napajet. Fro vypocet energie 
ulozene v akumulatoru mobilniho telefonu zvolme kapacitu K = 2 Ah a napeti U = 4,5 V. 
Potom pro ulozenou energii plati E = KU = 32,4 k.I. Kapacity baterii chytrych hodinek se 
pohybuji kolem K = 300mAh s napetim U = 3,8V, takze uchovavaji energii E = 4,1 k.J. 
Nakonec se podivejme na velmi prakticky pristroj, ktery se jiz vedci snazi napajet lidskym 
telem - kardiostimulator, jehoz prumerny prikon je 30pWT 

Turbm ka v krevnfm recisti 

Abychom mohli ziskavat energii z krevniho obehu, potrebujeme do nej nainstalovat nejaky 
system turbin. Z praxe vime, ze s rozmery turbiny roste i jeji ucinnost. Proto zvolme co nejmene 
turbin lined vedle srdce, kde jsou tepny a zily nejsirsi. Otazkou je z jake strany. Jestli do mista, 
kde krev ze srdce proudi (aorta a vystup plicniho obehu), ci naopak tarn, odkud krev proudi 
do srdce (liorni a dolin' duta zila, vstup plicniho obehu). K teto otazce se vratime pozdeji. 

Prace W, potrebna na j)rotlaceni objemu V, je soucinem sily a drahy, na ktere dana sila 
pusobi. Je-li na pocatku teto drahy tlak v kapalinc p\ a na konci tlak p 2 , vysledna sila bude 

F = S (pi - p 2 ) , 

kde S je obsah plochy prurezu trubice. Proudi-li kapalina rychlosti u, pro vykon bude platit 

P = Fv = Q {pi - p 2 ) , 

kde Q je objemovy prutok kapaliny v trubici. U zdraveho cloveka zhruba odpovida hodnote * 7 8 
Q = 0,071-tep" 1 = 0,084 hs” 1 . 

Tlak na zacatku turbiny p\ bude priblizne roven diastolickemu tlaku, ktery pro zdraveho 
cloveka uvazujme OOmmllg. Pokud budeme tok krve takto brzdit, bude sc srdce pravdepodob- 
ne snazit tuto ztratu dohnat tim, ze zvysi svuj vykon. Problemem je, ze tato zmena je velmi 
individualni (v zavislosti na zdravi daneho cloveka) a nelze ji prilis dobre predpovidat. Navic 
by tim dochazelo k nadmernemu zatezovani srdce. Pokud bychom vsak pocitali pouze male 
rozdily tlaku (pri malem odporu turbinky), mohli bychom tuto skutecnost zanedbat. Pro vetsi 
rozdily vsak muzeme vzdy udelat dolni odhad predpokladem pi = konst. Tlak p 2 umime nasta- 
vit zmenou odporu turbinky. Hledaine tedy co nejmensi diastolicky tlak, kdy jeste nedochazi 
k zadnym dalsim zdravotnim potizim. 

G https: //opentextbc. ca/anatomyandphysiology/chapter/24-6-energy-and-heat-balance/ 

7 https://link, springer, com/chapter/10.1007*/,2F978-3-642-50209-5_ 11 

8 https://hypertextbook.com/facts/2001/VitaliyShchupak.shtml 
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Takto zpusobena porucha se odborne nazyva ’’Isolated diastolic hypotension” 9 a nastava 
pri poklesu tlaku pod GOrnmHg, kdy clovek zacne pocitovat unavu a muze mit zavrate. Ty 
jsou zpiisobeny spatnym prokrvenim mozku, cemuz lze predcjit tim, ze dame turbinku na 
mista vtoku krve do srdce. Jedine nebezpeci tedy spociva v nadmernem zatezovani srdce a 
v problemech s tim spojenych, jako je napriklad zvysene opotrebovavani srdce. Ty se u lidi 
postizenych touto vadou projevuji az ve vyssim veku. 

Jaka bude ucinnost takove turbinky? Nami uvazovane zily maji prumer kolem 1 az 2 cm, 
coz jsou pro vyrobu prijatelne rozmery. Ucinnost vsak pri takovych rozmerech bude hral velkou 
roli, jelikoz s klesajicimi rozmery a vykonem vyrazne klesa. Pozorovanim zavislosti ucinnosti 
malych turbin na jejich vykonu 10 muzeme optimisticky odhadnout ucinnost vodni turbiny 
danych rozmeru na 10 %. Jelikoz ma krev priblizne stejnou hustotu jako voda (1060kg-m“ 3 ), 
nemusime tento odhad z duvodu pouziti jine kapaliny menit. 

Dosazenim hodnot do vztahu vysc dostavame 

P = r)Q (poo - poo) = 0,02 W . 

Nami stanovenou baterii telefonu bychorn timto zpusobem nabili za asi tri tydny. 

Pokud bychorn postupne zvysovali odpor turbinky (treba v ramci mesicu az let), mohlo 
by si srdce na vyssi potrebny tlak zvyknout a posilit se. Kazdopadne by tak dochazelo k jeho 
rychlejsirnu opotrebeni. 

Nektera zdravotni rizika jiz byla zminena. Co jsme vsak doposud nezminili je problem s niz- 
kou ucinnosti turbiny. Zbylych 90 % (ucinnosti jsme uvazovali jako 10 %) se totiz ve vysledku 
premeni na teplo. Toto teplo je vsak velini ryclile odvadeno krvi a pokud by nevznikalo tesne 
u sten zil, ktere by mohlo poskozovat, nemelo by byt nebezpecne. V posledni fade je treba 
zminit inoznou poruchu turbin, pri ktere by mohlo dojit k ucpani zil a nasledne vysoce prav- 
depodobne smrti. Dalsim velmi vysokym rizikem je vznik trombu kviili turbulencim, jez vedou 
ke krevnim srazeninam a nasledne srdecni zastave. V praxi bylo dosazeno 800 pW elektriny po- 
moci podobneho konceptu * 11 , coz je fadove mensi hodnota nez nami vypocitana. To muze byt 
zpusobeno spatnym odhadem ucinnosti nebo tim, ze jsme vykon turbiny tlacili az na samotnou 
hranici dlouhodobe prezitelnosti cloveka. Tu vsak lze redukovat snizenim odporu turbinky, cimz 
nastavime nizsi tlak, ktery pro srdce nebude pfilisnou zatezi. 

Mnohem slibnejsi koncept nez turbiny vyuziva nanotrubicek 12 a v budoucnu by jim bylo 
mozne napajet ruzne senzory uvnitr krevniho reciste, napriklad pro diabetiky nebo nanoboty. 

5/7a fuku 

Na podobnem konceptu by mohla fungovat turbinka umistena v hrtanu, kde vsak neproudi 
krev, ale vzduch. Potom muzeme pouzit vyse zmineny vzorec, upraveny pro danou situaci. 
Tedy vykon foukaci turbinky bude 


P — VQ (Pfuk Patm) • 

9 https: //www.uab. edu/mix/stori es/diastolic-blood-pressure-how-low- is-too-low 
10 https://ac.els-cdn.com/S187661021735124X/l-s2.0-S187661021735124X-main.pdf?_tid= 
b471df a5- 178d-4e57-92ee-a334313celc4&acdnat=1534083170_d5a9d5ddb0cd4230c8e48af 14375ae92 

11 http: //blogs .disco vermagazine. com/80beats/2011/05/17/tiny-turbine-inside-arteries-could- 
power-pacemakers-and-cause-blood-clots 

12 https://newsroom.wiley.com/press-release/angewandte-chemie-international-edition/ 
how-draw-electricity-bloodstream-one-dimension 
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Teraz to v rychlosti skontrolujeme so stredoskolskym riesenim. Rozlozime tiazovu silu 
v smere rovnobcznom a kolmoni so sikmou rovinou. Vidime, ze na teleso posobi len sinu- 
sova zlozka gravitacneho zrychlcnia 


a = g sin a , 

z coho je draha prejdena z pokoja za cas t pri tomto konstantnom zrychleni 


Aby sine to ale molili porovnat s nasim vysledkom vypocitanym cez Lagrangeove rov- 
nice, musime to trochu upravit. Nason zvolenou zovseobecnenou suradnicou bola :r-ova 
suradnica tclcsa. Chceme teda vediet iiosun telesa v x-ovom smere. Ten ziskame tak, ze 
ziskanu drahu prenasobime kosinusom uhlu ot. Teleso sa zaroveh pohybuje v opacnom 
smere ako je smer osy x, preto bude tento clen zaporny. Nazaver musime este pripocitat 
x-ovu siiradnicu vychodzieho bodu. Vychodzi bod je charakterizovany vyskou ( y-novou 
suradnicou) h. Musi preto platif 

h = xo tg a . 

Z toho jednoducho vyjadrime pociatocnu suradnicu a;o a dosadime do stredoskolsky zis- 
kaneho tvaru rovnice, cim dostaneme rovnaky vyraz, ako pri vypocte skrz Lagrangeove 
rovnice. 

Jakub Jambrich 
j akubj Qfykos.cz 
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3. Z hradiska toho, ako zvolit zovseobecnene suradnice, je tento problem triviality. Hmotny 
bod totizto kopiruje svojou trajektoriou linearnu funkciu zo sklonom a. Prevod medzi 
kartezskymi a zovsebecnenymi suradnicami bude preto 

Xm = x , 

Dm — X tg OL . 


Co nam dava analogickym sposobom ako v prvyeh dvocli ulohach vztah pre kinetickii 
encrgiu 

T — ^mx 2 (l + tg 2 a) . 

A takisto aj pre potencialnu energiu 


V = mgx tg a . 


Lagrangeova funkcia bude potom samozrejme 

L = ^mx 2 ( l + tg 2 a) — mgx tga. 

Napriek tomu, ze sine to v seriali explicitne nepreriesili na prfklade, bolo v nom odvodene, 
ako maju Lagrangeove rovnice vyzeraf, co nam dava navod, ako icli zostavit, ked’ inaine 
Lagrangeovu funkciu. V nasom pripade najprv parcialne derivujeme L podia x a vyslcdok 
nasledne podia casu. Od toho este odcitame derivaciu L podia x a cele to polozime rovne 
nule. Lagrangeova rovnica (bude samozrejme len jedna) ma teda tvar 

rnx (1 + tg 2 a) + mg tg a = 0. 


Rovnicu upravime 


x 


9 l g a 

1 + tg 2 a 


—g sin a cos a . 


Tento vztah staci dvakrat preintegrovat, cim dostanem 


x = s "‘ QrCOSa + C\t + C 2 , 


kde Ci ma vyznam pociatocnej rychlosti a C 2 ma vyznam pociatocnej pololiy. Ak po- 
ciatocna rychlost bola nulova a pociatocna vyska nad zemou h, pricom sme pociatok 
suradnicovej sustavy polozili do bodu, kde sa pretina naklonena rovina s podlozkou, po¬ 
tom v kartezskych suradniciach pre t — 0 z rovnice pre y mame 


h = x tg a , 
h 

tg a ' 


Kompletne riesenie tejto ulohy je 
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Ucinnost turbinky rj bude snizena diky tomu, ze ji proudi vzducli misto krve. Na druhou stranu 
bude zvysena tim, ze bude irioct byt vetsi. Vzhledem k tomu, jak hrube odhady zde delame, 11111 - 
zemc ji jiovazovat za stejnou jako v pripade krve. Objem vydechovaneho vzduchu a prislusny cas, 
ktery vydech trva, lze snadno dohledat. 13 Pro normalni nadech muzeme odhadnout V = 0,51. 
Jelikoz se nadechujeme pfiblizne sestnactkrat za minutu, turbinkou proteka Q = 0,13 dm 3 -s'” 1 
vzduchu. Pri dychani nase plice obvykle vytvareji pretlak 2 - 3kPa vuci okolnimu prostredi. 
Maximalm dosazitelny pretlak je pfiblizne 14 15 kPa. Odhadneme, ze pretlak 3 kPa oproti stan- 
dardnim 2 - 3kPa by nebyl prilis obtizny na udychani. Pak dostavame vykon P = 40 mW. 
Pro maximalm pouzitelny pretlak 12kPa (tedy navic oproti normalnimu pretlaku, nutnemu 
k dychani), mame P = IGOmW. Pomoci systemu potrubi a klapek bychom mohli dosahnout 
toho, ze by vzduch turbinou proudil stale stejnym smerem. Tim bychom docilili dvakrat vetsiho 
vykonu. 

Pro nizsi pretlak uvedme pro srovnani, ze nabiti mobilniho telefonu by trvalo temef 230 h, 
chytre hodinky bychom nabili za 30 h. Spotfebu kardiostimulatoru pokryvame pfiblizne tisic- 
krat. 


Odpadnf teplo 

Teplo, jez vyzafuji nase tela, ma z velke casti formu infracerveneho zafeni. Pro generovani elek- 
tricke energie uvnitf fotovoltaickych clanku je vsak nutne dosahnout urcite minimalni energie 
dopadajiciho zafeni. Typicka solarni kfemikova buhka vyzaduje alespoh 1,1 eV, takze funkc- 
nost je shora omezena vlnovou delkou dopadajiciho zafeni pfiblizne 1200 nm? 5 Zdola je pak 
omezena 10 hodnotou 300 11 m, viz graf 3. 

Uvazujmc absolutne cerneho cloveka, jehoz teplota je 37 °C. Cclkova intenzita zafeni bude 


■-*/ 


f 1200 nm 

I = n I dl (A), 

' 300 nm 


kde d/(A) je spektralni intenzita zafeni, kterou muzeme spocitat z Planckova vyzafovaciho za- 
kona. Nas vsak zajima elcktricky vykon na jednotku plochy, cili musime zapocitat jeste ucinnost 
solarnich panelu. Potom dostaneme 


Ps 


/• 1200 

./300 nn 


r,(A)d/(A). 


Numerickou integraci nam vyslo Ps — 45mW*m -2 . Stanovme plochu lidskeho tela 17 na 2nr. 
Vykon ziskany pomoci fotovoltaickych clanku z lidskeho tela tedy cim P = 90 mW, cimz muze¬ 
me nas telefon nabit za dobu pfiblizne 4 dm. Prakticky by bylo mozne clanky implementovat 
do obleceni, cimz bychom ziskali mnohem vice energie ze slunecniho zafeni nez z cloveka samot- 
neho. Dalsim nedostatkem fotovoltaickych panelu, jez nebyl zminen, je, ze pfi malych tocich 
energie nemusi vubec dojit k sepnuti nabijeci casti a ziskana energie je tak nulova. 


13 https://wwu.mada.org.il/en/about/engineer/challenge/respiratory-system 
14 http://flutopedia.com/refs_bpress.htm 

15 https://sciencing.com/effect-wavelength-photovoltaic-cells-6957.html 

1(i https: //www. researchgate.net/figure/Relative-quantum-efficiency-versus-wavelength-at-several- 

discrete-temperatures-for-a_figl_283165648 

17 https://www.calculator.net/body-surface-area-calculator.html?csex=m&bodyweight=85&bodyweightunit= 
kilogram&bodyheightf eet=&bodyheightinch=icbodyheight=180&x=96&y=23 
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Obr. 3: Ucinnost fotovoltaickych clanku v zavislosti na vlnove deice dopadajiciho zareni. 


i i-1-1-1-1-1-1-r 



Namerene hodnoty + 

Prolozena kfivka 

j _i_i_i_i_i_i_ ,j> m 


Vetrny pohon 

Obecne se udava, ze clovek za den vyprodukuje okolo 1,4 dm'* plynu v podobe vetrii , 18 jejichz 
slozeni je 19 

. Dusik: 20 - 90 % 

• Vodik: 0 - 50 % 

• Oxid uhlicity: 10 — 30 % 

• Kyslik: 0 - 10 % 

• Methan:0-10% 

Z tolio jsou horlave vodik a methan. Produkce methanu neni v lidskem tele nijak velika, 
avsak spalenim i jen maleho mnozstvi rnuzeme ziskat hodne energie, presneji 810kJ*mol~ l ?° 
Co se vodiku tyka, jeho spalenim lze ziskat az 120 k.)-g _1 neboli 240 k.J-inol -1 . Ten by take bylo 
mozne pouzit v palivovem clanku, coz je velmi perspektivni technologic s vyssi ucinnosti nez 
obycejne spalovani. 

Za normalnich podminek je molarni objem 22,4dm 3 -mol _1 . Za den tak dokazeme vyrobit 
energii maximalne 

W p = 12,6 kJ. 

Na vyuziti teto energie bychom potrebovali nejaky drub spalovaciho motoru, jejichz ucinnost 
se vetsinou pohybuje v intervalu 10 - 50%? 1 Uvazme tedy rj — 10 %, jelikoz i pri velmi dobrern 

1K https://www. telegraph. co.uk/men/the-filter/qi/10305094/QI-gas-facts-how-much-gas-does-the- 
average-human-produce.html 

19 https: //www. thoughtco. com/chemical-composition-of-f arts-608409 
20 https://www.wou.edu/las/physci/GS361/Energy_From_Fossil_Fuels.htm 
21 https://cs.wikipedia.org/wiki/Spalovaci_motor 
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suradnic pomocou zovseobecnenych suradnic 28 

X m = (/ + Z) COS if , 

Vm = {l + z) sin </?, 

Z M — Z . 

Casovym dcrivovanim jednotlivych suradnic najdemc vztahy prc jednotlive zlozky rych- 
losti v zovseobecnenych suradniciach. Nasledne urcime vztah prc kincticku energiu T 

T = (i 2 + (l -f z) 2 (p 2 ) -h t^Mz 2 . 

Nasledne urcime potencial. V tointo pripade sa jedna o potencialnu energiu zavazia (gu- 
licka ma pocas celeho pohybu rovnaku potencialnu energiu, nemusime ju preto uvazovat) 


V = Mgz . 


Lagrangeova funkcia pre tento problem bude potom vyzerat 

L = (i 2 + (/ + z) 2 + ^Afi 2 — Mgz . 

2. Znova si najprv urcime zovseobecnene suradnice. Kedze ulohu riesime ako rovinny pro¬ 
blem a vieme, ze trajektoria bodu M je zadana ako expiicitna funkcia (konkretne y = ax 2 , 
kde a je nejaka konstanta, aby nam sedeli jednotky), budu jeho obe suradnice vyjadre- 
ne pomocou jedinej suradnice x. Co sa tyka zovseobecnenych suradnic pre druhy bod 
spravajuci sa ako kyvadlo, mozeme pouzit suradnice, ktore sme si ukazali v seriali. Mu- 
sime vsak ale pamatat, ze jeho kartezske suradnice budu nielen suradnice kyvajuceho sa 
kyvadla, ale je nutne k nim pripocitat aj suradnice zavesu kyvadla o hrnotnosti M. 
Vztah medzi kartezskymi a zovseobecnenymi suradnicami je teda 

xm = x , 
y M = ax 2 , 
x in = x + l sin (p , 

JJm = ax 2 — / cos ip . 

Casovym dcrivovanim a dosadenim do vztahu pre kineticku energiu dostaneme 
T = i (M + m) x 2 (l + 4a 2 x 2 ) + \m ( cos *P + 2axsin <p) xip + / 2 ^ J ) • 
Potencial bude potom sucet potencialnych energii oboch telies, teda 
V = Mg ax 2 + mg (ax 2 — / cos (pj . 

Lagrangeova funkce bude jednoducho 

L = i (M -f m) x 2 (l 4- 4a 2 x 2 ) 4- (2/ (cos ip 4- 2axsin <p) x<p 4- l 2 p~) 

— Mgax 2 — rug (ax 2 — l cos . 

28 Tento ako aj nasledujuce kroky nebudem zvlasf komentovat, nakol’ko sa jedna o mechanicke pocitanie. 
Mudem len uvadzat, co robime a ako vyzera vysledok. Prcco to robimc (keby nickto ncvcdcl), je uvcdcnc 
v seriali. 
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Uloha III.S ... zobecnena 10 bodu; prumer 6,33; resilo 18 studentu 

1. Mcjmc vodorovnou dcsku, ve ktcrc jc mala dirka. Pres tuto dirku je provleccny provazek 
o deice l, na jchoz spodnim konci je zavcscno zavazi o hmotnosti M. Toto zavazi lzc 
povazovat za hmotny bod. Na druhem konci provazku na rovne desce je druhy hmotny 
bod (kulicka) o hmotnosti m. Provazek mczi nimi je napnuty diky zavazi o hmotnosti 
M. Cclou soustavu drzime v klidu tak, ze cast provazku pod dcskou je ve svislcm smeru. 
Pote druhemu hmotnemu bodu, kulicce, iidclimc rychlost v ve vodorovnem smeru kolmcm 
na provazek ve chvili, kdy soustavu uvolnime. V tomto pnkladu neuvazujte zadne treni. 
Zvolte vhodne souradnice a scstavte Lagrangcovu funkci pro tuto soustavu. 

2. Mejme zeleznou tyc ohnutou do tvaru paraboly tak, ze pokud v kartezske soustavc pusobi 
tihove zrychlem v zapornem smeru osy y, pak tyc ma stejny tvar jako funkce y = x 2 . 
Po tyci sc niuzc volnc pohybovat hmotny bod o hmotnosti M, ke kteremu jc pevnou 
nehmotnou tyckou o deice l pfipcvncno zavazi o hmotnosti m. Taktojsmc vytvorili kyvadlo 
se zavcscrn klouzajicim podel ohnutc tycc. Konstrukce dovolujc pohyb cole soustavy pouzc 
v rovine paraboly. Urcete vhodne zobecnene souradnice a najdete Lagrangcovu funkci teto 
soustavy. 

3. Mcjmc primku nakloncnou pod uhlcm a vzhlcdcm k vodorovnc rovine, po ktcrc sc po- 
hybujc bez treni hmotny bod o hmotnosti rn. Najdete vhodne zobecnene souradnice teto 
soustavy a sestavte Lagrangeovu funkci. Pote sestavte i Lagrangeovy rovnice, dvakrat 
je zintcgrujte, a tak najdete reseni. Zkontrolujte si, zda vase reseni vychazi stejne, jako 
reseni, ktere byste ziskali stredoskolskou metodou vypoctu. Pri integraci nezapomehte na 
integracni konstanty a vysvetlete jcjich vyznam. Jake budou jejich hodnoty, pokud se bod 
spusti z klidu z vysky h? 


1. Nazaciatok je dolezite urcit!, korko stupnov volnosti ma dana uloha, teda kolko zovseobec- 
nenych suradnic budeme potrebovat na jej popis. Vsimneme si, ze gulicka, pohybujuca sa 
po doske, ma dva stupne volnosti. Moze rotovat okolo dierky, cez ktoru je prevlecena nit, 
zaroven sa moze pohybovat v smere priarno d’alej alebo blizsie od dierky, nakolko spagat 
sa cez dierku tiez pohybuje bez trenia. 

Potom je tu zavazie, ktore tiez mozme reprezentovat ako hmotny bod. Ten ma ale len jeden 
stupen volnosti. Moze sa hybat hore a dole, inak povedane, ak zvolime osy x , y leziace 
v rovine, tak sa zavazie moze hybat iba pozdlz osi 2 (merane siiicrom nahor). Dalej ale 
mamc v tejto ulohe este jednu viizbu, tou je povrazok, ktory spaja hmotne body, a jeho 
dlzka je nemenna. Suradnica r, ktora udava vzdialenost gulicky od dierky, a suradnica 2 , 
ktora udava vysku zavazia, budu zviazane pomocou vztahu 

l — r — z . 

Na konstruovanie Lagrangeovej funkeie nam teda ostavaju uz len dve nezavisle suradnice, 
a to z a suradnica vyjadrujuca uhol, ktory zviera spagat medzi dierkou a gulickou s osou x 
(volba osi x je lubovolna vd’aka symetrii problemu). 

Teraz nasleduje druhy krok konstrukeie Lagrangeovych rovnic, vyjadrenie kartezskych 
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konceptu spalovaciho motoru nesmime zapoimnat na to, ze spalujeme pouze male mnozstvi 
plynu. Ucinnost spalovacich motoru totiz roste se zvysujici se teplotou. Takto by bylo mozne 
ziskat energii pribliznc W = 1,3 kJ. To staci na nabiti priblizne 4% baterie nami uvazovaneho 
telefonu. 

Seebeckuv efekt 

K tomuto jevu 22 dochazf pri spojeni dvou kovu ci polovodicu s ruznymi hodnotami Seebeckova 
koelicientu 22 a s rozdilnou teplotou. Mezi obema castmi vznikne napeti a tedy i elektricky 
proud. Napeti Ize jednoduse spocitat jako 

U = Ac* AT, 


kdc Aa je rozdil Secbeckovych koeficientu pouzitych materialu a AT je rozdil tcplot. 

Chceme vybrat dva materialy s co nejvetsim rozdilein Seebeckova koeficientu. V pripade 
polovodicu je nejvetsi rozdil mezi selenem a Pb 15 Ge 37 Se 58 , jez cini = 2890pV-K _l . Pri 
pokojove teplote 21 °C je teplotni rozdil vuci teplote lidskeho tela A'T = 16 K, coz odpovida 
napeti priblizne U = 46 mV. Jak je vysvetleno v clanku 23 o seebeckove koeficientu, nelze z polo¬ 
vodicu sestavit fungujici dratovy system a proto se pouzivaji spise pro opaeny Peltieriiv clanek 24 . 
Pro Seebeckuv je nutne pouzit kovy. Zvolme antimon a bismut, kde Aa = 119pV*K _1 , tedy 
U = 1,9 mV. 

Podivejme se na nejake realne koncepty. Napriklad v roce 2014 byl publikovan clanek 25 podle 
ktcreho bylo dosazeno vykonu na plochu priblizne l,5mW rn -2 . Jednim metrem ctvcrecnim 
bychom pokryli zhruba polovinu lidskeho tela. V livodu bylo zmineno, ze na nabiti mobilniho 
telefonu potrebujeme 32,4 kJ cnergie, coz by v tomto pripade trvalo zhruba 250 dm. Pro napajeni 
kardiostimulatoru by nam stacila plocha 200cm 2 . 

Tato metoda je velmi bezpecna. Co se tyce zatizeni organizmu, lidske telo je v zime pri- 
praveno zvysit svuj tepelny vykon z priblizne 100 W na zhruba 150 W a ochlazovani v radech 
miliwattu, dokonce ani wattu, mu nezpusobi zadne obtize, takze omrzlin se bat nemusime. Pro¬ 
blem by samozfejme nastal uz pri trochu vetsi zime, jelikoz kovove desky nelze povazovat za 
prilis dobre izolujici obleceni. 

Pro napajeni soucasnych mobilnich telefonu se vsechny tyto koncepty ukazaly velmi neprak- 
ticke. Pokud vsak uvazime podobnou elektroniku s mensim prikonem, napriklad energie ziskana 
z vetrneho pohonu ci krevni turbinky by pokryla polovinu elektricke spotreby chytrych hodinek. 
Navic vsechny tyto moznosti byly dostatecne na napajeni kardiostimulatoru, coz jc v dnesni 
dobe asi hlavni duvod zkoumani techto pristroju. Jako nejpraktictejsi se jevi vyuziti Seebeckova 
efektu ci vetrneho pohonu, jelikoz tato zarizeni by bylo mozne kdykoli sundat. Navic nemaji 
vedlejsi efekty, jako je napriklad mozna smrt v pripade turbinky v krevnim recisti. 

Jan Strelecek 
streldaOfykos.cz 

22 http: //vlab. amrita. edu/?sub=3tbrch=194&sim=351A;cnt = l 

23 https://www.electronics-cooling.com/2006/ll/the-seebeck-coefficient/ 

24 https://www.sciencedirect.com/topics/chemistry/peltier-effect 
25 http://iopscience.iop.org/article/10.1088/0964-1726/23/10/105002 
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Uloha III.E ... indexovany kondenzator 12 bodu; prumer 8,97; resilo 31 studentu 

Poston vam pfiscl clcktrolyticky kondcnzator a rczistor. Zmeftc kapacitu kondcnzatoru a odpor 
rezistoru, ncmcftc jc vsak pri'mo. Soucin kapacity kondcnzatoru a odporu r czistoru jc pribliz- 
neRCtt 20 s. 

Varovam: Elcktrolyticky kondenzator in a kladnou a zapornou elcktrodu, pri zapojem opacne 
ho muzetc znicit. Maximalm' dovolcnc napeti jc 10 V. 

Dodo prave v praktitku meHl rczonanci. 


Odpor rezistoru 

Najskor si zistime odpor rezistoru. Mamc ale zakazane merat ho priamo, nezostava nam teda 
nic ine, nez pouzit Ohmov zakon. Pripojime rezistor ku zdroju, odmeriame napatie na jeho 
koncoch a prud, ktory rrim preteka, a dosadime do znameho vztahu. 

Pri napati treba pripojit multimeter k tym dvom bodom obvodu, medzi ktory mi chceme 
vediet napatie, a teda paralelne k rezistoru. Kedze bateria uvadza 4,5 V nastavime na multimetri 
rozsah 20 V (pretoze 2 V by bolo zjavne malo a nic medzi tym nemame k dispozicii). 

Pri merani prudu zas treba, aby tento prud multimctrom priamo prctekal, a teda ho zapo- 
jime seriovo. Ak sa bojite, ze prepalite multimetru poistky, tak ho najprv nastavte na najvacsi 
rozsah prudu a postupne uberajte (zvysujte presnost), dokym sa da. 

My sme pri merani pouzili vzorku 5 rezistorov, pricom sa dalo cakat, ze mozno medzi nimi 
budu nejake drobne rozdiely, ale ukazalo sa, ze pri vsetkych rezistoroch sme dostali napatie 
U = 4,84 V a prud / = 96 pA. Co z toho vyplyva? Statisticka chyba odporov rezistorov vo 
vyrobe je zanedbatelna v porovnani s nepresnostou merania. Ta je v pripade nasho multimetra 
±( 1% + 2digit) pri merani prudu a ±(0.5% ± 2digit) pri merani napatia. Cize A/ = 3pA a 
A U = 0,04 V. Pri deleni dvoch velicin sa scitaju stvorce ich relativnych chyb. Takze odpor 
rezistoru je 



Trochu teorie 


Kapacita kondcnzatoru C udava pomer medzi nabojom Q v horn ulozenym a napatim U na 
jeho svorkach, teda C = Q/U. Podobne, odpor rezistoru R vyjadruje pomer napatia a prudu, 
ktory mm preteka, R = U/I. Co je to prud? Usporiadany pohyb nabitych castic. Jeho velkost 
udava, aky naboj presiel vodicom za jednotku casu, cize I = Q/t. 

Co sa teda stane, ked’ spojime konce nabiteho kondenzatoru rezistorom? Na svorkach re¬ 
zistoru bude rovnake napatie ako na kondenzatore a teda mm zacne tiect prud vel’kosti 



Q_ 
RC ' 


Tento prud tvoria clektrony bcziacc zo zaporne nabitcj elektrody kondcnzatora na tii druhu 
(resp. hypoteticke kladnc naboje bcziace presne opacnym smorom). Cize naboj a s mm teda aj 
napatie na kondenzatore bude postupne klesat. Tym padom l)ude klesat aj prud, cize dostavame 
diferencialnu rovnicu 

_dQ _ Q_ 
dt RC ’ 
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Diskusia 

Pruzky na rezistore maju postupne farby: zelena, hneda, cierna, cervena, hneda. Na internete 
sa da najst vela roznych moznosti ako icli dekodovat 27 Zistujeme, ze vyrobca udava, ze rezistor 
by mal mat odpor R 1} = (51,0 ± 0,5) kll. Tento interval je cely vnutri intervalu, ktory sme urcili 
nasim meranim, cizc hodnotu odporu sme urcili spravne. 

Vysledky kapacity ziskane meranim prudu a meranim napatia sa dost vyrazne lisia. Moze 
nas napadnut, ci to nemoze byt sposobene inou pouzitou metodou. Tak sme skusili aplikovat 
metodu merania napatia aj na meranie prudov. Prud klesol na polovicu (z 80 pA na 40 pA) 
v priemere za 18,2 s, co zodpoveda kapacite C — 520 pF (cize takej istej ako sme pri merani 
prudu urcili povodne), takze je problem niekde inde. 

Hodnota kapacity, ktoru udava vyrobca je C = 470 pF. Ta spada do intervalu chyby pri 
merani napatia, avsak nie pri merani prudu. To nam teda dava napovedu, ktore z merani je 
vlastne zle. 

A kde konkretne je ta chyba? Vratme sa spat k odvodeniu vztahu ( 1 ). Tam sme si povedali, 
ze R jc odpor, ktory stoji medzi svorkami rezistoru. Pri merani sme predpokladali, ze tento 
odpor je R = (50,4 ± 1,6) kO, co vsak nebola tak uplnc pravda. Seriovo s nasim rezistorom 
bol totiz zapojeny multimeter! Bol zapojeny ako ampermeter, takze jeho vnutorny odpor bol 
velmi maly (preto sme namerali skoro spravnu hodnotu), ale stale bol dost vel’ky na to, aby 
nezanedbatehie ovplyvnil meranie. 

Ked’ sa ale pozrieme na meranie napatia, vidime, ze aj tarn mohol multimeter ovplyvnit 
meranie. Tentokrat presne oimcnyin sposobom. Kedze bol zapojeny ako voltmeter paralelne 
s nasim rezistorom, skutocny odpor, cez ktory sa kondenzator vybijal bol mensi, nez sme po- 
citali, cize by sme mali pre zmenu dostat hodnotu kapacity nizsiu ako realna hodnota. Ale 
ostatne chyby v tomto pripade prevazili nad nedokonalostou voltmetru. 

S pouzitim d’alsieho multimetru by sme mohli odinerat vnutorny odpor v jednotlivych pripa- 
docli a pouzit ho na korekciu vysledkov. Do hry vsak vstupuju aj prechodove odpory v miestach 
zapojeni jednotlivych suciastick, ci vnutorny odpor samotneho kondcnzatora. Urcit presne ka¬ 
pacitu kondcnzatora tymito metodami je teda pomerne komj)likovane, najlepsie je totizto pouzit 
striedavy prud a zmerat induktanciu kondcnzatora, alebo urcit priamo rezonancnu frekvenciu 
kmitov prislusneho RC zapojenia. 

Zaver 

Hodnotu odporu rezistoru sme urcili ako R = (50,4 ± 1,6) kfJ, co sa zhoduje s hodnotou zadanou 
vyrobcom R v = (51,0 ±0,5) kfi. Meranim vyvoja prudu v case sme dostali kapacitu C = 
= (520 ± 20) pF. Meranim vyvoja napatia v case sme dostali kapacitu C = (480 ± 20) pF. 
Vyrobca udava C = 470 pF, co znamena, ze meranim napatia sme dostali spravny vysledok, 
zatial co meranim prudu sme dostali vyrazne chybny vysledok, dovodom bol nenulovy vnutorny 
odpor ampcrmctra. 

Simon Pajger 
legolasQfykos.cz 


27 https: //uww.digikey. com/en/resources/convers ion-calculators/conversion-calculator-resist or - color- 
code-4-band 


23 







Fyzikalni korcspondencm seminar MFF UK 


rocnik XXXII cislo 4/7 



4.5 i- 1 - 1 -i- 1 - 

+. 

4 - " 

3.5 - 

2.5 - - 

2 - 1 - 1 - 1 - 1 - 1 _ 

0 10 20 30 40 50 GO 

t 


Obr. 7: Linearizovana zavislost 


chyba smernice je J. s = 0,004. To spolu s relativnou chybou odporu dava uc = 0,03. Po 
vynasobeni kapacitou, dostaneme C = (520 ± 20) pF. 

Meranie napatia 

Tentokrat treba zapojit multimeter ku kondenzatoru paralelne. Rozhodli sme sa nabijat, dokym 
nebolo napatie 4,83 V, co trvalo asi sekundu. 

Dalo by sa samozrejme postupovat uplne rovnako ako minule a v pevne stanovenych casoch 
merat hodnotu napatia. Pre zmenu to ale spravime naopak. A sice, zmeriame casovy rozdiel 
medzi dvoma pevne stanovenymi hodnotami napatia. Tic sme si stanovili ako U\ = 4 V a U-i = 
= 2 V. Dolezity je hlavne ich pomer. Dosademm do ( 1 ) zistujeme, ze ak cas, ktory uplynie 
medzi tym, kym klesne napatie z U\ na U 2 , oznacime A/, mozeme C vyjadrit ako 


Spravili sme 5 inerani (opat jedno pre kazdy kondenzator) a urcili priemer a smerodajnu 
odcliylku At = (16,9 ± 0,4) s. 

Po dosadem do vztahu vyssie dostavame C = (480 ± 20) pF. 
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kde sme len za prud z predchadzajuceho vztahu dosadili zmenu naboja za cas. Minus je tain 
preto, lebo naboj 11 a kondenzatore tymto procesom klesa. 

Tato rovnica sa da riosit jednoducho separaciou premcnnych, ktorej vyslcdkom jc 

Q(t) = Q{ 0)e-7fc. 

Lcnzc naboj merat nevieme (minimalne bezne multimeter taku moznost neponuka). Zato 
vieme merat prud a napatie, cize mame dve moznosti: bud vydelit obe strany kapacitou, alebo 
obe strany zderivovat podia casu. Tym dostaneme casovy vyvoj napatia 11 a kondenzatore, resp. 
prudu tecuceho vodicom 

U(l) = U( 0)e"7fc , 

I{t) = /(0)e - ^ . (1) 

Vidirne, ze oba vztahy su analogicke, takze je vlastne jedno, ktory z nich pouzijeme. Ukazeme 
si samozrejme oba. 

Zostava otazka, ako vlastne zapojit obvod. Schemy zapojenia sa nachadzaju na obrazkoch 
4 a 5. 


" c 



Obr. 4: Zapojenie pri merani prudu 


Meranie prudu 

V tornto pripade rovno za kondenzator zapojime multimeter. Z predchadzajuceho merania uz 
vieme aky prud mame ocakavat. Multimeter teda momentalne ukazuje, ako rychlo sa konden¬ 
zator nabija. Tento prud tiez exponencialne klesa, takze je dolezite urcit si hodnotu prudu, pri 
ktorej nabijanie ukoncime. Nie je az tak dolezite aku, ale hlavne musi byt pri kazdom merani 
rovnaka. My sme sa rozhodoli, ze nabijanie ukoncime pri 3pA. Takto trvalo nabijanie nieco cez 
2 minuty. 

Ked’ donabijame, treba spravit sucasne dve veci - odpojit baterku (staci jeden kabel) a 
spustit stopky. Potom uz len sledujeme stopky a kazdych 5 sekund pozrieme na multimeter a 
poznacime si hodnotu prudu. Meranie sme vykonali pre pat kondenzatorov tak, ze sme tento 
postup zopakovali pre kazdy z nich, ale ako vidno z tabulky 2 , statisticka chyba je opat rozumne 
mala. 

Takze sme zmerali vyvoj prudu v case. Otazka znie, ako z toho urcit kapacitu. Najjednodu- 
chsi sposob je pouzit nejaky software, ktory nam nase data fitne exponencialou ? 0 

26 RStudio, gnuplot, . . . 
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R 


Obr. 5: Zapojenie pri merani napatia 
Tab. 2: Namerane casove priebehy prudu 


t 

s 



I 

pA 



i 

5 

78 

78 

80 

78 

79 

78,6 

10 

65 

65 

65 

64 

64 

64,6 

15 

54 

53 

54 

53 

54 

53,6 

20 

44 

44 

44 

44 

45 

44,2 

25 

36 

37 

36 

37 

37 

36,6 

30 

30 

30 

30 

30 

31 

30,2 

35 

25 

25 

25 

25 

25 

25,0 

40 

20 

21 

21 

21 

21 

20,8 

45 

17 

17 

17 

17 

17 

17,0 

50 

14 

14 

14 

15 

15 

14,4 

55 

11 

12 

12 

12 

12 

11,8 

60 

9 

10 

10 

10 

10 

9,8 


Okrem toho, ze sa mozcmc na grafe 6 pokochat, ako vsetky body lczia na dancj kriv- 
ke, sa dozvieme aj to, aky ma tato krivka prcdpis. V nasom pripade dostavame f(x) = 
= 94,3 • exp(—0,037 8x). Co keel porovname s (1) zistujeme, ze 


t 

lie 


-0,0378, 


RC = 26,5 s, 


C = 525 pF . 


Zaroveii si mozeine vsimnut, ze /(()) = 94,3 pA, co je o par mikroainperov inenej, nez priid, 
ktory sine namerali v prvej casti. Nahoda? Samozrejrne, ze nie. K tejto skutocnosti sa na zaver 
vratime. 

Ak ale nechceme fitovat namerane data pomocou nejakeho „black boxu“, mozeine to spravit 
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Obr. 6: Zavislost prudu v obvode na case pri vybijani kondenzatora 



aj sami pomocou linoarizacic. Ked zlogaritmujeme obe strany (1) dostavame, ze 

In/(f) = In 7(0) - j^ , 

cize ak si nakreslime graf zavislosti logaritmu ]>rudu (kedze do argumentu logaritmu by malo 
ist bezrozmerne cislo, budeme logaritmovat I/I n , kde sme si zvolili I n = 1 pA) na case, mal by 
byt linearny, so smernicou — ( RC )~ l , vid graf 7. 

Smernicu urcime z krajnych bodov 


A/(Q 

At 


-0,0379 s -1 . 


Z toho vyplyva, ze C = 524 pF, co je vlastne to iste, ako pri pouziti „black boxu“. Presnejsie 
by ale bolo vykonaf linearnu regresiu. 

Dalej by bolo vhodne urcif chybu. Stvorec relativnej chyby kapacity je sucet stvorcov re- 
lativnych chyb odporu (tu sme urcili v prvej casti) a relativnej chyby, s ktorou sme urcili 
smernicu. 

Ked si za oba prudy, ktore pouzivame na urcenie smernice dosadime / + u/, kde pod ui 
myslime absolutnu chybu daneho prudu, po chvili hrania sa s logaritmami a limitami sa da 
dojst k zaveru, ze absolutna cliyba smernice je 


Uk = 


u t 


( uh , 

Ui h ) ’ 


cize sucet relativnych cliyb prudov deleno casovy rozdiel. Ako absolutnu chybu prudu opaf 
pouzijeme dva krat najmensi dielik plus percento nameranej hodnoty a dostavame, ze relativna 


21 


















